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ABSTRACT 


A  method  is  presented  for  computing  stresses  in  the  vicinity  of  a  long 
cylindrical  cavity  in  an  infinite,  elastic,  isotropic,  homogenous  medium 
when  the  cavity  is  enveloped  by  a  plane  stress  wave  traveling  in  a  direction 
perpendicular  to  its  axis.  Hoop  stresses  around  the  cavity  boundary  were 
computed  for  the  passage  of  a  wave  of  dilation  and  a  wave  of  pure  shear; 
the  stresses  in  the  medium  away  from  the  boundary  were  computed  for  a 
wave  of  dilation.  Initially  the  stress  behind  the  incident  wave  front  was 
considered  to  be  constant.  The  effect  of  a  stress  decay  behind  the  front 
was  then  computed  by  using  the  Duhamel  integral  for  the  case  of  an 
incident  wave  of  dilation. 

The  method  of  solution  of  the  problem  involves  superposition  of  the 
stress  field  of  an  incoming  plane  step  wave  and  a  stress  field  corresponding 
to  waves  which  diverge  from  a  line  source. 
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I.  DESCRIPTION  OF  THE  PROBLEM 

The  purpose  of  this  investigation  is  to  determine  the  stresses  in 
the  region  of  a  cylindrical  cavity  in  an  infinite,  elastic,  isotropic,  homo¬ 
geneous  medium  when  the  cavity  is  enveloped  by  a  plane  stress  wave  traveling 
in  a  direction  perpendicular  to  its  axis.  Stresses  due  to  incoming  dilata- 
tional  waves  and  to  incoming  shear  waves  were  considered.  Initially  the 
stress  behind  t)ie  incident  wave  front  is  considered  to  be  constant,  but  some 
consideration  is  given  to  the  effect  of  pressure  decay  behind  a  shock  front 
through  the  use  of  the  Duhamel  integral.  The  problem  is  one  of  plane  strain; 
that  is,  each  section  perpendicular  to  the  axis  of  the  cavity  is  in  the  same 
state  of  stress,  there  is  no  displacement  in  the  direction  of  the  cavity 
axis,  and  only  three  stress  components  must  be  computed.  Some  consideration 
is  given  to  the  variation  of  the  parameters  which  define  the  properties  of 
the  medium. 

The  practical  problem  corresponding  to  the  dynamic  elasticity  pro¬ 
blem  described  above  concerns  the  effect  of  stress  waves  initiated  by  nuclear 
explosions  on  underground  protective  installations.  Although  any  practical 
problem  of  this  type  cannot  exactly  reproduce  the  conditions  of  the  theoreti¬ 
cal  problem,  it  may  be  reasoned  that  under  certain  conditions,  particularly 
under  those  for  which  such  protective  construction  is  likely  to  be  built, 
fairly  good  agreement  may  be  expected.  For  example,  rock  is  neither  iso¬ 
tropic  nor  homogeneous,  but  a  protective  structure  would  be  constructed  at 
a  site  where  the  rock  is  as  sound  and  unstriated  as  possible.  The  assump¬ 
tion  of  an  infinite  medium  is  more  reasonable  than  it  might  at  first  appear. 
For,  in  this  case  a  protective  structure  would  surely  be  several  diameters 
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telov  the  ground  surface,  and  since  the  maximum  hoop  stress  is  reached  at 
about  four  to  five  transit  times  after  the  wave  reaches  the  opening,  reflec¬ 
tions  from  the  surface  will  not  interfere  with  the  peak  stress.  Probably 
the  most  severe  restriction  is  that  of  linear  elasticity-  However,  many 
rocks  approach  this  property  in  the  lower  range  of  their  stress-strain 
diagrams.  The  actual  opening  is  not  likely  to  be  a  perfect  cylinder;  never¬ 
theless,  current  mining  procedures  can  produce  an  opening  which  approximates 
a  cylindrical  shape  quite  well. 

The  importance  of  this  problem  has  been  recognized  by  several 
earlier  investigators.  In  Ref.  (l)  a  study  was  made  of  the  stresses  in  the 
medium  using  a  "high  frequency"  approximation  similar  to  the  approach  used 
in  geometrical  optics.  This  method  does  not  allow  the  solution  to  be 
carried  out  for  sufficient  time  to  obtain  the  maximum  stresses.  Reference 
(2)  considered  the  stresses  around  a  hole  In  a  thin  plate  when  enveloped 
by  a  plane  stress  wave  of  harmonic  variation  in  time  and  space  behind  the 
wave  front,  i.e.  the  plane  stress  problem  corresponding  to  the  plane  strain 
problem  treated  herein.  Reference  (5)  presented  a  solution  to  the  problem 
investigated  in  the  present  report  obtained  by  using  an  integral  transform 
in  time  of  the  displacement  potentials.  The  resulting  equations  were  then 
solved  in  terms  of  Hankel  functions.  Considerable  difficulty  is  encoun¬ 
tered  in  this  method  in  the  computation  of  the  inverse  transform.  In  the 
present  work  an  entirely  different  method  of  solving  the  partial  differential 
equations,  which  has  several  advantages  over  the  previous  method,  is  presented. 
The  results  of  the  two  analyses  will  be  compared  in  Chapter  VI.  In  Ref.  (4) 
the  method  of  Ref,  (3)  was  used  to  solve  the  same  problem  again  and  was 
extended  to  include  stresses  away  from  the  boundary.  Also  possible  failure 
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mechanisms  "based  solely  on  the  elastic  stress  state  were  discussed.  Results 
were  presented  in  the  form  of  contour  maps  of  principal  stresses  at  each 
time  interval.  The  angle  of  minimum  principal  stresses  and  angles  of 
possible  failure  planes  are  also  presented.  In  Ref.  (5)  the  solution  to 
the  same  problem  was  formally  carried  out  by  a  method  quite  similar  to  that 
used  in  Refs.  (3)  and  (4).  The  theoretical  derivation  of  the  equations 
was  reported  in  Ref.  (5),  but  numerical  results  had  not  yet  been  obtained. 

The  problem  has  also  been  considered  experimentally.  In  Ref.  (6) 
a  combination  of  photoelasticity  and  grid  analysis  was  used  to  determine 
the  stress  distribution  on  the  boundary  of  holes  in  a  plate  of  low  modulus 
material.  The  shape  of  the  holes  and  duration  of  stress  wave  were  varied. 

The  method  results  in  a  stress  field  of  both  dilatational  and  shear  waves 
with  a  rather  long  rise  time  at  the  front  of  the  wave.  Reference  (7) 
contains  the  results  of  a  large  number  of  experiments  made  by  detonating 
an  explosive  charge  at  the  surface  of  blocks  of  brittle  material  which  con- 
tained  openings  of  various  shapes.  The  medium,  shape  of  opening  and  distance 
from  a  free  surface  to  the  opening  were  varied.  These  experiments  are 
useful  in  predicting  the  type  of  failure  to  be  expected. 
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II.  METHOD  OP  ANALYSIS 

The  "basic  method  of  solution  of  the  problem  involves  superposition 
of  the  stress  field  of  an  incoming  plane  step  wave  and  the  stress  field 
corresponding  to  waves  which  diverge  from  a  line  source.  The  diverging 
waves  are  chosen  so  that  the  radial  and  shear  stresses  due  to  the  two  waves 
combine  to  produce  a  traction-free  cylindrical  boundary  centered  on  the  line 
source.  This  cylindrical  boundary  is  then  the  boundary  of  the  cavity. 

A  solution  to  the  scalar  wave  equations  corresponding  to  a 
cylindrical  wave  diverging  from  a  line  source  may  be  found  in  Refs.  (8)  and 
(9)  •  This  solution  can  be  applied  to  the  dynamic  elasticity  problem  by  con¬ 
sidering  the  displacement  potentials  to  be  represented  by  Fourier  series 
expanded  in  6,  with  the  coefficients  written  as  functions  of  radial  dis¬ 
tance  and  time.  These  series  are  substituted  into  the  wave  equations, 
resulting  in  hyperbolic  partial  differential  equations  in  two  independent 
variables  which  must  be  satisfied  by  the  coefficients  of  each  term  of  the 
series.  The  general  solution  of  these  differential  equations  is  then 
written  in  integral  form;  therefore,  the  boundary  conditions  result  in 
integral  equations. 

In  order  to  satisfy  the  boundary  conditions,  the  sum  of  the 
stresses  due  to  the  diverging  wave  and  the  incident  stresses  expanded  in 
Fourier  series  were  set  equal  to  zero  at  a  fixed  radius  from  the  line 
source.  In  general  this  results  in  two  linear  integral  equations,  one  from 
the  radial  Stress  condition  and  one  from  the  shear  stress  condition,  which 
must  be  solved  simultaneously  for  two  unknown  functions  appearing  in  the 
integrands  of  the  integrals  mentioned  above.  The  integral  equations  were 
solved  numerically  in  order  to  determine  these  unknown  functions . 
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A  program  for  solving  the  simultaneous  Integral  equations  and  com¬ 
puting  the  various  stresses  at  the  boundary  and  in  the  medium  avay  from  the 
boundary  was  written  in  Fortran  language  for  solution  on  a  CDC  1604  digital 
computer.  This  program  will  be  described  in  Appendix  C. 

In  order  to  check  the  results  obtained  from  the  computer,  the 
same  basic  equations  were  solved  by  another  method.  Each  of  the  Fourier 
coefficients  of  the  displacement  potentials  for  the  diverging  wave  and  the 
incident  stresses  was  expanded  in  a  Taylor  series  about  zero  time.  By 
considering  only  the  first  few  terms  of  the  series,  it  was  feasible  to 
carry  out  the  computations  by  hand.  The  results  thUB  obtained  are  only 
valid  when  the  variable  in  the  Taylor  expansions,  time,  iB  small.  For  this 
reason  the  method  is  called  the  short-time  solution.  This  additional  solu¬ 
tion  is  useful  as  a  check  on  the  computer  solution  at  early  time  and 
also  gives  the  Initial  values  of  the  unknown  functions  in  the  integral 
equations,  which  are  needed  in  the  machine  solution.  As  a  further  check 
on  the  computer  solution,  the  static  stresses  on  the  boundaries  and  in 
the  medium  were  compared  with  those  obtained  from  the  machine  solution  at 
very  long  time. 

This  general  method  of  solution  was  applied  to  the  problems  of 
an  incident  wave  of  dilatation  and  an  incident  wave  of  pure  shear.  Fur¬ 
ther  details  of  the  representation  of  the  incoming  stress  field  and  of 
the  stress  field  diverging  from  a  line  source  are  described  in  Sections 
5.1  to  3.4.  The  equations  resulting  from  the  combination  of  these  two 
waves  to  obtain  the  boundary  equations  are  derived  in  a  form  suitable  for 
use  in  a  congmter  in  Sections  4.1  and  4.2. 
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The  diverging  wave  which  has  been  determined  at  the  boundary,  of 
course,  continues  to  travel  outward  and  becomes  the  reflected  and  diffracted 
waves.  To  find  the  stress  in  the  medium  away  from  the  boundary  the  same 
superposition  procedure  may  be  used,  except  that  the  effect  of  change  of 
radial  distance  of  the  diverging  wave  must  be  considered.  The  computer 
solution  for  stress  in  the  medium  due  to  a  plane  wave  of  dilatation  is 
presented  in  Section  4.3. 
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HI.  DERIVATION  OP  EQUATIONS 


-  1  Fourier  Analysis  of  Incident  Wave  of  Dilatation 

First  consider  a  plane  step  wave  of  dilatation  traveling  in  the 
negative  x  direction,  with  a  stress  ex  1-  the  direction  of  wave  propaga- 
non  It  may  he  shown  that  the  stresses  m  the  y  and  z  direction  are 


ffy  =  TI^y  ax 

az  =  7I7TT  Jx 

where  v  is  Poisson* c  ratio  For  convenience  v  will  be  used  as  defined  by 


v 


"Ti^T  * 


3-1 


In  the  following  work  tensile  stresses  will  be  considered 

positive 

Consider  the  stresses  cn  a  circle  of  radius  a(in  cylindrical 
coordinates)  being  enveloped  by  the  wave  front  described  above,  with 
rv  --  c,  and  with  no  cavity  o.n  tr.e  ncdium.  Referring  to  Fig-  >  .1  the  radial 
stress,  arr,  shear  stress  ?  0 ,  and  he  op  stress  behind  the  front  nay  be 
written  using  tie  genera",  rule.,  ci  trarst'er. .aticn  of  stress-, 


drr  -  <s  Uo-."  •?  ♦  i  - n"  ’-■>> 

(l-\) 

c  „  - -c  * — ~ —  -*ir‘ 
rB 

o.a  =  i  (-in"  9  +  i  cos'  9) 

UO 


Sine?  o  and  <y,,.  are  even  m  9  while  cr  is.  edd  in  9,  the 
rr  ed  rB 

Fourier  series  corresponding  to  those  expressions  can  be  w-nttem 
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a  (t)  r~> 

a  =  — 75 —  +  )  a  (t) 

rr  2  /  .  rr  1 

n=l 

oy 

ar0  =  X  bn(t)  sin  n6 


cos  n0 


<*  (t) 


°QB  ~  2  + 


I  dn(t) 


cos  n0 


The  coefficients  in  these  series  may  he  evaluated  from  the 
following  integrals  where  0^  defines  the  extent  to  which  the  wave  has 
enveloped  the  circular  boundary  ac  shown  in  Fig.  3.2,  and  therefore  is  a 


function  of  time 


2  r1 

an(t)  =  -  J  crrr  (0)  cos  n0  d0 
o 

8  r 1 

bn(t)  =  —  /  0^(0)  sin  n0  30 

o 

2  r1 

dn(t)  =-  J  aQe  (0)  cos  nfi  dfi 


If  time  is  measured  from  the  instant  that  the  wave  touches  the 
boundary,  the  plane  diletational  wave  with  a  velocity  of  propagation 
will  travel  in  the  time  t 


from  which 


Cjt  =  a  (1  *  cos  0^) 


ex  =  coo  U~tC1/a) 


0  <  tCL/a  <  2 


0L  =  n 


tC^/a  >  2 
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Carrying  out  the  integrations  in  Eqs.  one  obtains  the  follow¬ 

ing  expressions  for  the  Fourier  coefficients  as  a  function  of  time: 

For  *  <  la,  C, 


a  (t) 
o  *  ' 


■  f  [h  +  si»  aeJ 


3.6 


rsin(2~n)0.  sin(2tn)0 , 

”„<*)  ■  -  S  ^-^L'Ta-nT--  -  “ISXr 


•] 


3.8 


d«(*) 


g 

2  ~  2x 


^(l+v)©1  +  sin  20l 


r  sin  n0,  /,  sin^2+nJ0,  sin(2-nJ0. 

„<*>  ■  f  If1"'  — 5-9 


~v  sin(2+n)0.  sin(2-n)01 


For  t  >  2a/C, 


ao(t)  g(l+v) 

2  ~  2 

&n(t)  =0,  n/2,  a2(t)  - 

bn(t)  -  0,  n  s*  2,  l2(t)  =  - 

do(t)  _  q(l+v) 


1.10 


dn(t)  =  0,  n  {  2, 


4,(0  =  - 


g(l  v) 


Equations  3*3  and  3  ,  combined  with  Eqs-  ; .t  to  3 .10  give  the 
time  dependence  of  the  free  field  incident  stresses  on  a  cylindrical  boundary 
with  radius  a.  The  variation  of  the  first  three  of  these  coefficients  with 
time  is  shown  in  Fig.  3 . 3  for  a  value  of  v  =  1/3  and  with  a  =  1  0.  The 
accuracy  with  which  the  first  three  terms  of  a  Fourier  series  represent  the 
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stress  variation  with  angle  around  the  opening  at  two  times  may  be  observed 
in  Fig.  3.4.  In  this  figure  the  variation  of  a  ,  a  and  aaa  with  angle 

is  shown  for  the  times  tc^/a  =  0.5  and  1.0.  Also  shown  on  this  figure  is 
the  actual  variation  of  these  stresses  for  comparison.  It  may  be  observed 
from  these  curves  that  the  representation  is  not  as  good  for  the  earlier 
time  as  the  later,  as  would  be  expected.  For  this  reason  the  method  pre¬ 
sented  is  not  applicable  for  early  stress  confutations  unless  a  larger 
number  of  terms  is  used  in  the  Fourier  series.  The  representation  for  the 
incoming  wave  is  exact,  however,  for  all  times  after  the  wave  has  enveloped 
the  boundary.  Also  it  may  be  seen  that  the  representation  is  much  better 

for  a  and  a  „  than  it  is  for  a00.  Since  the  diverging  wave  form  is 
rr  rt/  y  y 

determined  by  the  incident  boundary  stresses  cr  and  a  „ ,  the  determination 

rr  tv 

of  the  diverging  wave  should  be  fairly  accurate. 

3.2  Fourier  Analysis  of  Incident  Shear  Wave 

Consider  the  case  of  a  plane  step  wave  of  pure  shear  with  magni¬ 
tude  of  stress  t  traveling  in  the  negative  x  direction.  Stresses  on  a 
cylindrical  boundary  behind  the  front  may  be  written, 

a  =  ?  sin  20 
rr 

=  ?  cos  20  _  3*44 

3qq  -  -  Y  sin  20 

where  positive  shear  is  defined  in  Fig.  55-  A  tilde  will  be  placed  over 
those  quantities  associated  with  an  incident  shear  wave  in  order  to 
distinguish  them  from  the  corresponding  quantities  associated  with  an  inci¬ 
dent  Wave  of  dilatation.  The  Fourier  series  corresponding  to  these  expres¬ 
sions  axe. 
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CTrr  =  bn(t)  sin  ne 

n=l 

a  (t) 

5r0  =  ~^2 —  +  an^  cos  li0 
n=l 

CO 

°ee  =  ^  5n(t)  sln  110 

n=l 


3.12 


The  Fourier  coefficients ,  found  in  the  same  manner  as  for  waves 


of  dilatation,  are 
for  t  <  2a/C0 


for  t  >  2a/C- 
— - 


b  (t)  =  - 
nv  « 

rsin(2-n)Sg 

sin(2+n)0g-i 

L  (2-n) 

(2+n)  J 

a  (t)  =  - 

n  it 

rsin(2-n)02 

sin(2+n)0o-| 

L  (2-n) 

+  (2+n)  J 

dn(t)  -  - 

?  rsin(2-n)02 

sin(2+n)02 

it  [  (2-n) 

(2+n) 

bn(t)  =  0, 

n  1  2, 

bg(t)  =  t 

Sn(t)  -  o. 

n  *  2, 

a2<  t)  =  r 

d  (t)  =  0, 
rr  ’ 

n  1  2, 

cU(t)  =  -  T 

% 


3.15 


3.14 


where  C0  is  the  velocity  of  propagation  of  a  shear  wave.  The  angle  ©2 
may  be  written  as  follows  by  referring  to  Fig.  5*.1 


for 

and 

for 


eo  =  cos'1  (1  -  tC2/a) 

0  <  tC2/a  <2.0 

0.,  -  n 

£l 

tC2/a  >  2.0  . 
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3 - 3  Derivation  of  Equations  for  Diverging  Wave  (incident  Dilatational  Wave) 
The  equilibrium  equations  may  be  expressed  in  the  following  form 
for  an  elastic  medium 

(X+p)  +  p  u  »  p  u  =s  0 

(X+p)  ^  +  p  -  p  v  =  0  3*15 

(X+p)  ^  +  pA  -  p  ii  =  0 


where  A  is  the  dilatation  given  by 


A 


5u  5v  5w 
=  5^  +  &  +  SI 


and  u,  v,  and  w  are  displacements  in  the  x,  y  and  z  directions.  It  is 
convenient  to  introduce  the  following  representation  for  the  displacement 
vector  u  in  terms  of  a  scalar  potential  q>  and  a  vector  potential  f 


u  =  grad  9  +  curl  f  . 


Ihen  grad  9  represents  the  irrotational  part  of  the  displacement  and 
curl  5  represents  a  displacement  corresponding  to  no  dilatation.  Observing 
from  the  physical  problem  that  there  is  no  variation  in  displacement 
with  z  along  the  cavity  axis  and  that  there  is,  therefore,  only  rotation 
about  the  z  axis,  one  may  write  the  following  simplified  expressions 


grad  9  =  ^  +  ^ 

curl  5  =  ^*  -g?<N 
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Therefore 


•setting  5Z  = 

The  displacements  u,  v,  and  w  may  "be  written  then 


w  =  0 


The  equations  of  motion  (Eqs.  3*15)  will  he  satisfied  if  the  potentials 
<p  and  v  satisfy  the  following  two  wave  equations, 


f  =  cp 

V2  ?  =  V> 


3.17 


where  the  velocities  of  the  dilatational  and  shear  waves  are  given  by 


X  +  2p 
P 


and 


5.18 


Ir.  polar  coordinates  the  displacement  components  of  Eqs.  5 .16  are  given  by 


u  -  I  _  &£ 

0  r  §0  or 


3.19 


where  u  is 
r 

tive  in  the 


taken  positive  in  the  direction  of  increasing  r  and  u^ 
direction  of  increasing  0 


posi- 
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In  Ref.  (8)  the  following  general  solution  to  equations  of  the 
type  represented  by  Eqs.  3-17  is  given.  First  the  potential  (displacement 
potentials  in  this  case),  tp  and  y,  are  taken  in  the  following  form: 


r11  f  (r,t)  cos  n9 
n 


V  = 


I 

n=1 


rngn(r,t)  sin  nfl  . 


3.20 


These  particular  forms  for  the  displacement  potentials  were  chosen  because 
they  satisfy  ^qp  =  0  and  V^y  =  0,  a  "homogeneous"  part  of  the  differential 
equations  if  fR(r,t)  and  gn(r,t)  axe  assumed  constant.  These  functions  of 
r  and  t  must  be  chosen  so  that  the  equations  with  qJ  and  V  on  the  right  are 
satisfied.  Since  the  displacements  u^  caused  by  the  incident  wave  of  dila¬ 
tation  are  even  in  0  while  the  displacements  u0  are  odd  in  the  0,  <p  is 
chosen  even  and  y  odd  so  that  Eqs.  3*19  will  result  in  displacements  of 
the  same  form  when  cp  and  y  are  substituted  into  them.  Therefore  this 
form  for  representing  the  diverging  wave  is  only  suitable  for  use  with 
an  incident  wave  of  dilatation. 

Substitution  of  Eqs  3  20  into  Eqs  3.17  gives  the  following 
differential  equations  to  be  satisfied  by  the  Fourier  coefficients, 

2  ,S  fn  2n+l  Sfn.  „ 
ci  +  —  SF5  =  Fn 


Co  (- 


V 


2n+l  ^aii. 


5.21 
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The  solution  to  these  equations  may  be  written  in  terms  of  the  solution 
for  n  =  0  as 

f  -(-■%-)  ?  and  g  =  (-^-)  g  3.22 

n  'r  or'  o  °n  vr  dr'  °o  •  ^ 


In  Section  B.l  of  Appendix  B  it  is  shown  that  Eqs.  3-22  do  indeed  satisfy 
Eqs  3.21.  The  next  step  is  to  find  tQ  and  g^.  This  may  be  done  in  the 
following  way.  Consider  the  first  of  Eqs.  3.21  with  n  =  0. 


a  a2f  .  3f 

<4  5?>  - 

This  is  recognized  as  the  wave  equation  in  cylindrical  coordinates  for 
the  case  of  symmetry  about  the  z  axis.  The  form  of  f  which  satisfies 
Eq  3  23  may  be  obtained  through  the  following  reasoning.  It  is  known 
that  the  form  of  the  displacement  potential  for  a  spherically  symmetrical 
point  source  wave  is 


The  displacement  potential  for  a  line  source  may  be  found  by  integrating 

the  displacement  potentials  for  all  point  sources  uniformly  distributed 

along  the  z  axis  from  z  =  -«  to  z  =  +«,  each  point  source  having  the  same 

variation  with  time.  Then  the  symmetrical  line  source  displacement 

potential  may  be  written 
+x 

ro  -  /  S  f  (t  -  c:>  d*  ■  5'24 

-OC  1 

This  same  form  of  displacement  potential  is  used  in  Ref.  9  in 
the  investigations  of  supersonic  fluid  flow  around  a  long  wing- 
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5.28 


17 


From  Eqs.  3  .28  and  Eqs.  J.20, 


CO 

n  / 1 

r  (F 

+OQ 

J)  /  F(t  - 

n=0 

0 

CO 

+oo 

rn  (i 
vr 

5?)  /  G(t  " 

n-0 

0 

Considerable  simplification  can  be  effected  in  the  solution  by- 
noting  that  Eqs.  J.29  can  be  written  in  the  following  form 

00  +03 


tp  =  f"  y  f  l,a(t  -  ~  cosh  u)  cosh  nu  dul 

kfco  «x"  J0  '  J 


cos  n 0 


r  oo  +oo 

^  =  '  X jf  (t  -  ~  cosh  u)  cosh  nu  dujsin  n B} 
n=0  C£  0  2 

where  the  superscripts  on  F  and  G  indicate  the  order  of  the  derivative. 

The  proof  of  the  conversion  from  Eqs.  J.29  to  3. JO  is  shown  in  Section  B.2 
of  Appendix  B. 

In  terras  of  displacements,  the  stresses  at  any  point  in  the 
medium  may  be  expressed  as  follows: 


a  >A  +  2p 
rr 


°r8  =  u  'r0 


aee  "  ^  +  f '  +  ur} 


bur  u  ?u0 
A  ^5F+T  +  er  7 
au0  t  dur  u0 
7r9  =  57"  *  7  59  V 


By  substituting  into  these  equations  the  displacements  given  by 
Eqs.  3 “19 j  one  obtains  the  stresses  in  terms  of  displacement  potentials. 
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O  =(x  +  2u)^£  +  ^i£+i^^  +  ?H.  _  §&  |£ 

rr  ^  2  r  §r  +  2  -n2  +  r  5r§9  2  08 

or  r  o#  r 

rS  ~  H  ordS  2  §0  -.2  2  -..2  r  dr  1  -5’-52 

r  or  r  do 

,  32(P  .  (X  +  2u)  5q>  .  (X  +  2u)  d2ep  2u  b%  ,  2p  cty 

aee  -  x  Z  2  +  ?  5r  +  “2  io2  “  T  5?§e  +  ^  5e  • 


Substitution  of  the  general  expressions  for  9  and  (Eqs.  3.50)  into 
Eqs,  3.32  gives  the  general  expressions  for  stress  in  the  diverging  wave. 
00  u, 

\  1  f/  i\n  p ^ ^  o 


cos  n0 


=  ^ J  Fn+2(q1)[(2u  cosh2u  +  X)  cosh  nu]  du 
n=0L°l  0 

-  £  Gn+«-^  ^  cinh  2u  sinh  nu]  duj  coo  n0 

00  n  U1 

=  ^  J'  Fn+2(q1)[d  oinh  2u  sinh  nu]  du 

n=C>  L°1  0 

-  Gn+c(TVo)  td  cosh  £u  cosh  nu]  duj  sin  n0 

C2  t> 

co 

=  X  Fn+2^ni^  ^?v  ■ 2pi  sinh2u^  °o3h  nu^ du 

n=0  L°1 

n  u2  "I 

+  ■{— f  Gn+t(rw)  tu  sinh  Cu  sinh  nu]  du  cos  n9 

r  n+^  -L  *  J 


2_>  L n  vA 

n=C>  °1  0 

(-l)n  f2  Gn+: 

■  c  n+2  J  G 

b.3  0 

co  U_ 

_  y  riii^  r 

L  n+2  J, 

n=0  “1 


C,  ‘•0 


where 


’b.  =  (*  “  <J"  cosh  u)  and  ^  =  (t  -  jr* 


cosh  u) 


3.36 
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The  integration  Units  on  Eqs.  3.24  and  therefore  on  Eqs.  3. 30 
extend  to  +»  since  there  are  point  sources  all  along  the  z  axis.  The  wave 
emanating  fron  each  point  source  is  considered  to  have  the  same  variation 
in  time.  Therefore,  at  any  finite  time  only  the  point  source  waves  out  to 
a  finite  value  of  z  have  reached  the  point  in  question  provided  they  all 
started  at  the  same  finite  time.  For  a  finite  time  the  integration  in  u 
need  only  extend  to  the  value  corresponding  to  this  value  of  z. 

By  substituting  Eqs.  3-30  into  Eqs.  3*19>  the  general  expressions 
for  displacements  may  also  be  found. 


pi 


V1  T  (-l)nr  f1  r.n+2,  wsinh  (l+n)u  sinh(l-n)un  .  , 

Ur  =  L  L  i  F  (V[  i?liny  +  . "^t^-]sinh  u  du 

n-0  L1 

*  f  -  4==  * 


3-37 


u 

a(-l)nr  /  „n+2,  Nrsinh(l+n)u  sinh(l-n)u-,  ,  , 

-  r^rj  F  (tii)[  Ta^r — 2^r-]  sinh  u  du 

n=0  4  0 

*  I*  — » 41”  -  • 


These  egressions  for  stress  and  displacement  in  a  wave  diverging 
from  a  line  source  are  only  suitable  for  elimination  of  the  boundary  stresses 
due  to  an  incident  wave  of  dilatation  since  the  variation  in  9  has  the  same 
fori,  in  the  two  cases. 


3  A  Derivation  of  Equations  for  the  Diverging  Wave  (incident  Shear  Wave) 

The  displacements  in  the  medium  caused  by  a  plane  incident  Bhear 
wave  with  a  constant  state  ox'  stress  behind  the  front  are  odd  in  8  in  the  radial 
direction  and  even  in  the  0  direction-  In  order  that  the  displacements 
caused  by  the  diverging  wave  (Eqs.  3*2-9)  might  have  this  same  variation 
<p  should  be  odd  and  $  even  in  8.  The  displacement  potentials  in  this  case 
are  then  the  same  as  for  an  incident  wave  of  dilatation  except  for  the 
variation  in  0,  and  may  be  written  as  follows  by  analogy  with  Eq.  3.29: 

00  +00 

$  -  Y,  r^r  5r^n  I  F(,1l)  du  sln  110 

n=0  ^ 

00  +00 

?  =  y  rn(i;  |p)  /  G(q2)  du  cos  n0  . 

n=0  v 

These  expressions  may  be  written  by  analogy  with  Eq.  5*50 


ffl  =  y  f  Fn(ri, )  cosh  nu  du]  sin  n0 

n4okn  J0  1  J 
00  +00 

5  -  I  AF  /  0%)  cosh  nu  duj 

n-n  Lco  u 


3.38 


cosh  nu  duj  cos  n0  , 

Substitution  of  Eqs.  3*38  into  the  general  expressions  for 
stress  in  terms  of  displacement  potentials  (Eq.  3*32)  gives, 


rr 


cosh^u  +  \)cosb  nu)  du 
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The  two  equations  resulting  from  these  boundary  conditions  are  to  be  solved 
for  the  functions  F11'^^)  and  G114^^), 

For  an  incident  wave  of  dilatation  the  boundary  equations  may 
be  found  by  combining  the  first  two  of  Eqs  J with  Eqs.  3.33  and  3-3*i 
and  applying  the  resulting  conditions  at  the  radius  r  =  a. 


0  =  an(t)cos  n8  n+2  J  Fn+2(q1)  [(2^  cosh2U4\)cosh  nu]du 

n=l  n=0  L°1  0 

J  Gn+2( T)g)[tJi  sinh  2u  sinh  nu]duj  cos  n0 


vn  pi 


r  ^ 
0 


n  U1 

0.  ^bn(t)  P kf  sinh  2u  sinh  nu]  du 


n=0  "1  0 


~~n42  J'  G^^iigHn  cosh  2u  cosh  nu]  du  sin 


“2  0 


In  a  similar  manner  the  boundary  equations  for  an  incident  shear 
wave  may  be  found  by  combining  the  first  two  of  Eqs.  3.12  with  the  first 
two  of  Eqs.  3 .39. 


0=  l  t>n(t)  sin  n 9  4  J'  Fn+“(i^M2u  cosh2U4>.)cosh  nu]  du 

n=0  n=l  C1  0 

,  ,n  p2 

4  ■' ~  ■  /  Gn+£:(ru)(|i  sinh  2u  sinh  ru]  du  oin  n£? 

n42  J  *2 


a  (t) 

°  =  -v+ 


u. 

n  pi 


^  an(t)  cos  n0  ^H+2  f  Pn+£‘(\)t"dSinh  2u  sinh  nu]du 


n  G-, 
n=0  1 


22 


The  two  equations  resulting  from  these  boundary  conditions  are  to  he  solved, 
for  the  functions  Kn+I_(  t^)  and  Gn+W(q2) , 

For  an  incident  wave  of  dilatation  the  boundary  equations  may 
he  found  hy  combining  the  first  two  of  Eqs,  3.3  with  Eqs.  3.33  and  3.34 
and  applying  the  resulting  conditions  at  the  radius  r  =  a. 


0  =  -T5 —  an(t)cos  n0  +  ^T  J  Fn+2(  tj^)  [(2p  cosh2u+\)cosh  nu]du 

n=l  n^6  Lci  0 

n  u2 

-  —g  J  Gn+2(Tig)[u  sinh  2u  sinh  nu]duj  cos  nfl 


\n  A. 


\n  o2 


n  ^ 

0  =  ^bn(t)  sin  n0  +  jsj  J'  sinh  2u  sinh  nu3 

n=l  n=0  L°1  0 

-  ^  J'  Gn+"(  tjg) [u  cosh  2u  cosh  nu]  duj  sin  n0 


"2  0 


In  a  similar  manner  the  boundary  equations  for  an  incident  shear 
wave  may  be  found  hy  combining  the  first  two  of  Eqs.  3*12  with  the  first 
two  of  Eqs.  3.39. 


0=  X  **n^  sin  n9  +  ^ J'  Fn+"(Ti^)[C:u  cosh2u+>.)cosh  nu]  du 
n=0  n=l  C1  0 

+  f  Gn+2(tjo)(u  sinh  hu  sinh  ru]  dul  sin  n© 

r  n+t"  J 

0 

a  (t)  “  n  U1  '  *42 

0  =  —^5 —  +  an(t)  cos  +  ^  j  ~"^2  f  F°+2(  [ -qsinh  2u  sinh  nu]du 

n=l  n=0L°l  ^ 

n  U2 

-  J'  Gn+2(r^)  [p  cosh  2u  cosh  nu]  duj  cos  n 9 


0  =  -^-  + 


(-l)n  r%  n+2 


G0  'b 


The  number  of  parameters  involved  in  the  solution  of  the  problem 
may  be  found  from  Eqs  J.4l.  The  last  two  terms  in  each  equation,  which 
give  the  stresses  in  the  outgoing  wave,  contain  only  X  and  p  as  parameters 
But  X  may  be  expressed  as 

.  2v 

x  =  T&  * 

By  making  this  substitution  for  X,  one  may  factor  p  out  of  both  terms, 
leaving  v  as  the  only  parameter  that  cannot  be  factored  out.  The  actual 
form  of  an(t),  bn(t),  and  dR(t)  may  be  seen  by  referring  to  Eqs.  3 -6-3 .10, 
where  it  is  found  that  the  only  parameter  within  each  term  is  v,  a,  the 

incident  stress,  being  a  common  factor  of  each  term.  Since  v  is  a  func¬ 

tion  cf  v  (Eq  3.1),  Poisson’s  ratio  is  also  the  only  parameter  which 
cannot  be  factored  out  of  these  terms.  In  the  actual  machine  computations 

all  terms  were  divided  by  p  and  the  ratio  set  equal  to  unity,  i.e.  p 

was  taken  as  the  unit  of  stress.  Then  a  single  parameter  appears  in  the 
solution,  v.  In  the  derivation  of  equations  which  follows,  X  and  p  are 
retained  in  the  equations  so  that  it  is  clear  to  what  extent  the  variation 
of  these  constants  alters the  solution  of  the  problem. 

All  the  constants  mentioned  thus  far  occur  in  the  equations  for 
hoop  stress  and  equations  for  the  stress  in  the  medium  in  precisely  the 
same  way  us  mentioned  above.  Therefore,  the  same  comments  may  be  made 
concerning  the  parameters  involved  In  the.se  computations. 

The  sane  observations  say  be  nade  concerning  Eqs.  5 .*+2  for  an 
incident  shear  wave,  except  that  t  is  a  common  factor  in  each  tern  for 
stress  due  to  the  incident  wave  (Eqs.  j  lb  and  $.1U)  and  v  does  not  appear 
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The  last  two  terns  of  each  equation  are  the  sane  as  before  and  again,  p 
may  be  factored  out  leaving  v  as  the  only  parameter  that  cannot  be 
factored  out.  In  the  computer  solution  't/p  was  set  equal  to  unity  in 
this  case. 

In  the  numerical  solution  of  the  problem,  the  fundamental  unit 
of  length  was  taken  as  a,  of  velocity  (C-,  for  incident  shear  wave),  and 
of  stress  p.  These  three  units  completely  determine  all  the  other  units 
required.  In  Eqs.  3 -^l  the  value  of  C,  is  unity  while  the  value  of  C„  is  a. 
Where  a  is  defined  to  be  the  ratio  q  .  and  are  nevertheless 
carried  through  the  derivation. 
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IV.  MODIFICATION  OF  EQUATIONS  FOR  COMPUTER  SOLUTION 

4 . 1  Equations  for  an  Incident  Wave  of  Dilatation 

Consider  the  numerical  solution  of  the  boundary  equations  given 
by  Eqs.  3-4l  for  an  incident  wave  of  dilatation.  The  first  terms  on  the 
right  involving  a^(t)  anil  bn(t)  represent  the  stresses  due  to  the  incident 
plane  vave .  These  terms  are  functions  of  time  and  may  be  evaluated  by 
computing  the  time-dependent  angle  0^  from  Eq.  J.5  and  substituting  this 
angle  into  the  appropriate  expression  for  the  Fourier  coefficients  given 
by  Eqs.  $.6  -  5.10. 

The  last  summations  in  Eqs.  5.41  represent  the  stresses  in  the 
diverging  vave  and  contain  integrals  which  must  be  evaluated  to  obtain  the 
stress  It  is  important  to  note  that  these  integrals  are  of  the  convolu¬ 
tion  type  and  therefore  the  integrations  must  be  performed  anew  for  each 
value  of  time  considered.  At  any  time  for  which  the  stress  in  the  diverg¬ 
ing  wave  is  to  be  evaluated  the  integration  can  be  performed  numerically 
provided  all  the  previous  values  of  Fn+“(iu)  and  Gn+“(ip,)  are  known.  It 
is  these  functions,  however,  which  must  be  found  from  the  boundary  condi¬ 
tions-  Thi  .  is  possible  since  Zq.  .5.41  constitutes  two  equations  with 
two  unknowns. 

In  the  numerical  integration  process  it  is  convenient  to  take 
equal  _teps  in  the  independent  variable-  However,  if  equal  steps  in  u^ 
and  u-  are  taken ,  unequal  steps,  in  time  and  unequal  steps  in  movement  of 
the  outgoing  wave  in  sp.  t-e  will  result.  To  avoid  this  -.c.rplication  the 
variable  of  integration  will  be  changed 
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Referring  to  Fig.  3.6  and  the  equation  following  Eq.  3. 85  one  may- 
write  the  limits  on  the  integrals,  u^  and  Ug  for  the  time  t, 

r  cosh  =  p1  =  r  +  tC^ 
r  cosh  Ug  =  p2  =  r  +  tCg. 

Solve  for  cosh  u^  and  cosh  Ug 

tCl 

cosh  u^  =  1  +  — p-  =  1  + 

tc2  .  4,1 

cosh  Ug  =  1  +  =  1  +  Cg. 

Therefore  Eq.  3«4l  may  he  written,  using  the  substitution 


cosh  u  =  1  +  £  , 


4.8 


with  the  limits  ^  on  the  integrals  involving  FIi+2(t>l)  and  ^  on  those 
involving  Gn+2(T^).  These  upper  limits  may  he  expressed  as  follows 
from  Eq.  4.1  for  the  particular  case  r  =  a. 


tc. 


5i- 


and 


fe¬ 


te. 


By  using  the  substitution  of  Eq.  4.2,  one  may  write  the  arguments  of 
Fn+S( r^)  and  Gn+2(T|g)  given  hy  Eq.  3.36 

ru  =  t  -  — -  (1  +  0  and  Tfe=t-?r-(l+£).  ^-3 

1  t1  ^  Ug 

In  perforaing  the  numerical  integration  it  is  convenient  to  have  the 
limits  on  all  the  integrals  the  same.  This  will  he  accomplished  hy  chang¬ 
ing  the  upper  limit  on  the  integrals  of  Gn+2(rjg)  to  agree  with  that  for 
j,n+2(  )  1>e>  both  Upper  Units  hecane  To  do  this  the  variable  of 
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integration  for  the  integrals  involving  Gn+2(  t^)  must  be  multiplied  by  the 
ratio  of  the  limits 


h 


-  a. 


4.4 


For  simplicity  the  subscript  on  the  upper  limits  will  be  dropped.  Then  the 
upper  limits  for  all  the  integrals  is 


tc. 


or,  for  the  particular  case  r  =  a, 
tCl 

7  =  — i 


4.5 


and  the  'variable  of  integration  f  for  the  integrals  involving  Gn+"( r^) 
will  to  replaced  by  at,  ■ 

In  the  remaining  explanation  of  the  numerical  solution  of  Eqs.J.41 
it  is  convenient  to  consider  a  particular  value  of  n  in  order  to  simplify 
the  equations  involved.  The  case  n  =  1  will  be  considered,  however,  the 
derivations  are  quite  similar  for  n  equal  to  any  other  value  except  for 
certain  simplifications  when  n  -  0 

Consider  the  integral  portions  of  Eqs.?  41  representing  stresses 
due  to  the  diverging  wave  with  n  -  1  (these  stresses  are  also  given  by 
'Ll..'  and  ’■  ?-) 


r^L.  f 1 

1  7.  I 

c. " 

L  J 
ir- 

77  I 


F'  ( rp  ) [ ( Ce  cor.h'u  *\)  cosh  u]du 
G  *(i).){u  si  nr.  lu  sinh  u]  dul  cos  0 
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h'ur.erical  integration  of  these  equations  is  hindered  by  the  fact 

that  one  term  of  each  stress  equation  has  an  integrable  singularity  at  £ 

equal  to  zero  It  is  also  necessary  to  assume  the  functions  and 

q- ’  (i^)  to  be  well  behaved,  however  this  assumption  is  justified  for 

snail  t  by  the  short  time  solution  which  will  be  discussed  in  Section  5*1* 

Also  there  is  an  integrable  singularity  in  F* • ’ ( at  the  value  of  the 

argument  it  =  +a/C^,  however  it  turns  out  that  this  does  not  affect  the 

stresses  obtained  from  the  numerical  integration.  Consider  as  an  example 

the  numerical  integration  of  the  first  term  on  the  right  of  Eq.  4,6.  The 

l/2 

method  of  integration  to  be  used  assumes  all  of  the  integrand  except  £  ' 

in  the  denominator  to  be  linear  between  discrete  values  of  £.  Weighting 

l/2 

factors  consistent  with  a  linear  variation  in  £  divided  by  £  are  com¬ 
puted  and  applied  to  the  ordinates  at  each  end  of  the  interval  in  order  to 
find  the  value  of  the  integral  in  the  interval.  The  derivation  of  the 
integration  formula  is  shown  In  Appendix  A  This  method  of  treating  inte¬ 
grals  of  the  type  encountered  here  is  similar  to  a  method  derived  in 
Ref  (10)  If  £  is  divided  into  equal  steps,  A£,  and  the  above  idea  applied, 
the  value  of  the  integral  between  m  A£,  and  (m  +  l)Af,  for  the  first  term 


on  th*'  right  of  Eq  -  6  may  be  written 
( r,+l  )A£ 

(F'  ( VJr!uV.  ['^r.mm+l)1/'-^]  ft(m+l)'/2-m',/S] 


F  l\) 


(h-t-lt.)  (v*-CQr  *  1 


it 


=  -  JSf 


*  <r  "C*U  ]  -  » 
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and  m  is  the  numbering  of  steps  in  At  starting  with  zero  at  £  =  0.  For 
brevity  the  following  symbols  are  used  for  the  weighting  factors  in  Eq.  4. 9. 


Am  -  2(l+m)[(m+l)'!'/‘~  -  j  [(m+l)''^2  -  m^*"] 

Bm  =  |  [(m+l)5y/2  -  m5/2]  -  2m[(m+l)'1‘^2  -  m1/2]  . 


4 .11 


It  is  shown  in  Appendix  A  that  these  weighting  factors  may  be  used  for  any 

of  the  integrals  involved  in  the  stress  equations  as  long  as  the  function 

to  be  integrated  is  well  approximated  by  a  linear  relation  m  £  divided 
l/2 

by  £  '  in  any  interval. 

It  would  only  be  necessary  to  use  the  special  integration  scheme 
for  the  first  few  steps  near  t  -  0;  however,  in  the  computer  solution  it  is 
convenient  to  continue  to  use  it  throughout  the  range  of  integration.  The  total 
integral  may  be  expressed  as  the  summation  of  all  the  integrals  in  the  interval 
A£  as  given  by  Eq.  *.9, 

The  following  abbreviations  for  quantities  in  Eqs.  u .6  to  t .8  are 

used: 
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Ri'  _  j" -t- I 

01  "  I-  VT  c\  J 

rG  =  f  WC(l^)(l-^Zg?|  a 
01  *-  V 2a  N !  1  tat/2 


The  stress  equations  given  by  Eqs.  4.6  -  4.3  may  be  written  as 
f\  1  Lews  for  any  time,  1  =  pAt ,  where  p  =  0,1,2,...,  by  using  the  summations  of 
Eq,.  4. Q  for  the  integrals  and  the  abbreviations  of  Eqs.  4.11  and  4.12 

r  P-1 


rr 


I  ^  [(f"  ■  (v>f-s(Bii)«v(F’ "  (,ii)  wAa] 

n=0 

r^1  r  -i" 

>  \ (G,"('lJ)  (r°  )  a  +(g’ "(O)  .(R*.)  ,b 

i  I  '  '  l2//r-r.  Alb  m  '  v  ‘2"p-m-l'  Al'm+1  mj 


m-0 


COS  6 

4.13 


rti 


P-1 

ir*-0 

p-l 

£  -tot  f(0"'!v»)t.,(^1)lrAn.W"h2))p.,.1<4)»AJ 


r.-<' 


r-i 


sin  0 


4.1; 


) 


r-i 


•'  ‘p  -r,<*01  »r.V(F’  '  '  (T‘l  *  Wl<  *DlW 


r.  -0 


cos  0 


-».i' 


In  order  ‘•hat  the  value  ♦’  E  -ay  te  "r.p  i+ti  tr.ee  ar.i  used  earn  tire  tne  inte¬ 
grations  are  performed,  eqaal  .steps  In  tire  are  .'hssen  and  steps  At,  are  used  as 
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given  by 

C1 

AC  =  At  —  . 

*  a 

The  use  of  Eas.  4,15  and  4.14  in  numerical  fern  permits  the  boundary- 
conditions  given  by  Eqs.  J.41  to  be  rewritten  as  follows  for  the  case  n  =  1. 

In  these  equations  the  unknown  values,  (F‘  1  ( rjj > >p  and  (G1 '  '(^))p>  axe  removed 
from  the  summations. 

0  -  (a^lUj,  <0'  <’fe>V8«4Ao] 


p-2 


I  '  ''\)W,-i(RL4.r<0' ■  <  VW 

n=0 

<o-<v.v5i>,]  Vi 


4.16 


o » (vt))p^  [(F,,,(^v4i1oV(a',;{’>2,v4vAcJ 

P-2 

*1  ■^[(p,',<’>iWil4w  •  (O'  -id J  (VVi> 

m=0 


4.17 


These  equations  are  to  be  applied  at  r  a. 

The  stress  caused  toy  the  incident  wave,  ie.,  a^(t)  and  b^ft),  is  a 
known  function  of  time  and  therefore  may  be  evaluated  at  times  p&t.  If  all  the 
previous  values  of  F!,(,q^)  and  <3 '  ' ' 1  q  . )  an  kn  vr,  t  r.  n  tf  .  \  ulu*-v  of  the 
summations  are  computed  except  for  the  terms  involving  (F',!(q.  ))  and  (G '  ' '  (  T|0  '  . 

1  p  p 

The  values  of  these  unknowns  may  then  be  found  from  the  resulting  equations. 

Since  the  first  value  of  time  that  can  be  considered  with  these  equations  is 
t  =  At,  it  is  necessary  to  find  the  values,  (F'!  iq^j}^  and  (G' '  Uq^))^,  by  some 
other  means.  Then  the  values  at  t  =  At  can  be  computed.  Knowing  the  values  of 
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T.he  functions  at  t  -  0  and  t  =  At,  those  at  t  =  2At  may  he  computed,  and 


:o>  on. 

Xr.e  short  tir.e  solution  ty  which  the  initial  values  of  F' 1 '  (np)  and 
d‘"  (''i-,)  were  found  is  discussed  in  Chapter  V,  where  the  values  are  given  as 


(F’"(n;L))0  =  -2cJa/[jt(\+2p)] 

(G'"(n-,))-  -  o. 

c.  u 


For  machine  solution  Eqs.  4.l6  and  4.17  may  he  written 

(F--*(ni))p(Rli)0A0  *  (G'"W)p(Rj1)0A0  =  - -  (sA1)p 

(M*)) 


(*''*<VWoV 


-  (s  ) 

s/ZT  B1  p 


4.18 


4.19 


where 


P-2 


"  2.  [tF"  '  (''l))p.r..ltEAl4+l  +  UTi3))p_n_jCRA1)m+1](Bm+Am+1) 

*  «f">\)4(EL4  *  “"'''.'“AV  Vi] 


r.-O 


4.20 


P-2 

"  L 


Or.’-,  uli  ’ to-  vaiie.-  F’  ur.d  G '  "  ( t  , )  are  known  ->p  *  ‘  a  given 

n.rVl  i_.  f-o'-irl*-  *•:  fir. i  ’h>.  h*  r-p  o  t  re  o  'it  thi.  tir.e.  Ir.  order  perform 
t.v.  sdrrvatitns  for  o  y  along  with  the  ether  summations  it  is  necessary  to  write 
4.1=  f:  r  the  res*  die  to  the  diverging  wave  in  tne  save  form  ’■hat.  Up.  t.lj 
and  -  were  writ*-.-.,  -=-.'ir.;r  trv  ‘err.  involving  (?",(ri1'l)p  and  (3*  **  (»»,)) 

out  of  the  sumr.at  j O.ue  the  unto,  -i  .  .  r-.  f-uni,  ‘he  ’or-  tun  le  aided  V 


is  added 


the  suorat ion .  Also  the  stress  due  to  tr.e  incident  wave,  (d^(t) )^, 
to  that  caused  by  the  diverging  wave. 

•*>  ■  *  (q"' h2>  ao] 


[  (*"  "  (Tl1))p.m.1(RD;L)n;+1  +  (G‘ ' 1  (n2) )p.m_i(BD]L)m+;L]  (b0+Ab+1) 


m=0 


4.21 


^[(F'"(V)0(i)p  +  (G,,,(n2))0(Ro1)  ]  B 


To  compute  the  incident  stresses  at  any  time  it  is  desirable  to  have 
the  variable,  t,  in  a^t),  b^Ct),  and  d^(t)  written  in  terms  of  t  =  pAt.  This 
may  be  done  by  dividing  the  time  required  for  the  wave  to  travel  one  radius 
into  N  steps  and  starting  time  when  the  plane  wave  first  reaches  the  side  of 
the  opening.  Then  the  ratio  p/lJ  is  equal  to  tC^/a  where  p  is  the  step  number 
corresponding  to  the  time  t  =  pAt,  until  the  wave  reaches  the  back  of  the 
opening.  Therefore,  Eq.  5*5  may  be  written 

0^  -  cos  1  (l  -  p/H)-  4.22 

Once  0^  is  known  Eqs.  3.6  through  5. 10  may  be  evaluated.  For  the 
case  n  ^  1  and  p  <  2N  these  equations  became 


*i(t)  =  if  [(5  +  v)  sin  di +  siri  H 


bq(t) 


C.TX. 


(1  -  v) 


tin  t* 


-  %  [<5v 


1  -  |  =ir' 

•  “f1  :  -Mi]  • 


4.23 


and  for  p  >  2H 


ax(t)  =  0,  bx(t)  =  0,  d1(t)  =  0. 


4.24 


The  case  n  =  0  is  considerably  simplified  by  the  absence  of  a  shear 
wave  throughout  the  computation.  There  is  only  one  equation  then  to  compute 
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for  the  easi=>  n  -  ^  are  as  follows: 


The  equat  i?:. 


T,T  I*  tv  r*  (a.(t)) 

: F1  ‘  (t,  '  1  (<  _)  A  +  (GU(tu,)  (R1;,)  a  =  -  --—■>■  -  (S*  ) 
v ‘1  -fv  Aa  o  o  v ‘2  7p  A-  o  :■  /rr  A2  f 


lfIV'W4>o  A0  ,  (G-(r,,))  (^)o  Au  = 


IMO) 


£  _  /oP 


*  (s^)T 


•**>,  - 1  [A>w4Wi  *  (‘"'VwA.ij  <vw 

nr*=0 

*  |yV<N  »:<»!;•'>(  ’  t»,V<-.h  ),(»*-!,] 

t-2 

«V,  -  I  [(FlV<''l))I-r-l<4)r..l  •  <sIV(\»r-r.-1(4)B.lJ  <VW 

r.-0 

+  [(fIV(\>VR£A  *  Bp_, 


4.26b 


*  (GIV(T,2))p(R°2 


tPo  Aj 


^  A  f^VWl'4',..  •  ‘^W'kW] 

r.  0 

[(FIV(rtl));fRFJr  V  'GiVO.J).(r:J.)rl  Br_, 


A*  i-  c;v4  J 

rg  _  _  ra^cu^c/G)1^^)^  . 


a* 4 1 
<£  ^2  -J 
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si.  -  -  1  a 

'-a  l*a£/2  -* 

,F  _  + 

“K  I-  cjVjsT 1^71  j 

!f]a 


4.26d 


a*c£  •/i# 


(FIV(ri,))-  -  +4C*a/[(>.+2ti)«],  (GIV(i|2))o  =  0 


4.26e 


f-t  f  <  2:: 


>  4.26f 


u^(t)  ~~  [(i-v)  '■>-  +  (i+v)  sir.  2o,  +  t-  (l-v)  sin  40.] 
t.(‘)  -  -  [*.  -  ~  sir.  4o  ] 

i.C)  [Cv-1)  +  :i  +  v)  si:,  .o,  +  4-  (v-1)  sin  46^  ] 


. r  i  is. 

M--v),4,  ty  )  -  -o(:-v)/2,  !>(♦.)  *  -o(l-v)/2  4.2cg 
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boundary  are  small  and  go  to  zero  at  one  transit  tirr.e.  Since  the  stresses  in 
the  outgoing  waves  are  determined  by  these  incident  stresses,  it  is  reasonable 
to  expect  that  they  should  also  be  small  and  die  off  quickly  after  one  transit 
time  similar  to  the  stress  for  n  =  1,  which  also  goes  to  zero  at  the  back  of 
the  opening.  The  final  stresses  ccmputed  by  using  the  first  three  modes  will 
be  presented  and  discussed  in  Chapter  VI. 


4.2  Equations  for  Incident  Shear  Wave 

The  numerical  solution  of  the  bcundary  equations  for  an  incident 
shear  wave  (Eqs.  3.42)  is  quite  similar  to  the  procedure  discussed  above  for 
an  incident  wave  of  dilatation.  It  may  be  seen  from  Eqs.  3-13  and  3.14  that 
the  incident  stress  on  the  cylindrical  boundary  dees  not  depend  on  Poisson's 
ratio.  It  may  be  observed,  however,  from  Eqs.  3-42  that  the  stresses  caused 
by  the  waves  which  diverge  from  a  line  source  do  depend  on  Poisson's  ratio 
since  these  equations  contain  both  >.  and  p. 

It  is  again  convenient  to  consider  the  solution  for  a  particular 
value  of  n  in  order  to  simplify  the  explanation;  the  value  n  =  1  will  be  used. 
Also  since  the  derivation  is  so  similar  to  the  case  just  considered,  the  dis¬ 
cussion  will  be  somewhat  abbreviated. 

Consider  only  the  integral  terms  m  Eqs .  3  42  which  reprerorst  the 
stresses  in  the  diverging  wave,  and  let  n  -  i  (also  given  by  ‘he  first  two  cf 

Eqs  3.59). 

’"1 

°rr  =  H  /  F"’(n1)[(2p  ccsh^u+>.)  cosh  u]  du 


1 

I 


G'  sinh  2u  sinh  u]  du  I  sin  t? 
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i 

°rd  "  f  "T  I  F‘  *  *  C*li)tp-  sinh  2u  sinh  u]  du 

^  -j 

J  G’  ''(iigJtji  cosh  3u  cosh  u]  duj  cos  0 


4.27 


c\  “o 


Also  the  hoop  stress  in  the  diverging  wave  may  he  written  by  setting  n  =  1  in 
the  last  of  Eas.  3*  39- 


00 


4  /  F1 1 '  (tj  )[2h  sinh2u  1  >-)  cosh  u]  du 
Jo  1 

u2 

-4  /  G' ''(t|0)[h  sinh  2u  sinh  u]  du  sin  0 
P*  Jr,  J 


4.28 


0'  ~o 


Again  the  substitution  for  cosh  u  in  terms  of  £  as  given  by  E4.  4.2  will  be  used, 
uni  the  upper  limits  on  the  F'  ’ '  (r)1 )  and  G' ' '  (r,^)  integrals  become 

~  +  C1  ~  tC2 
f  —  and  £,  «  —  . 

S1  a  a 

In  this  case  it  is  convenient  to  have  the  upper  limit  on  the  integrals  involving 
F' 1 '  ( T|  )  agree  with  that  on  the  integrals  involving  G"'(il2).  In  order  to 
change  the  upper  limit  to  on  the  F' ' '  (t(,  )  integrals,  the  variable  of  Integra - 
•  i  r,-  r.us*  he  cult  ini  iei  by  the  ratio  of  the  old  to  the  new  limits, 

L-fi.i., 

t/cr° 

For  simplicity  the  subscript  will  be  dropped  on  the  upper  limit;  then  the  limit 
for  f he  integrals  becomes 

x  r<t 

o'  - 

*  =  1T  * 

Tne  above  epuatic.s  for'  stress  r.ay  then  be  written  as  follows: 
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The  jane  numerical  integration  pro  -eiare  tlvt  was  use!  earlier  in 
♦his  chapter  uni  ier*vei  in  Appendix  A  nay  te  applied.  Evrationa  4.2?  nay  be 
written  -ulv*  w-ustaticns  for  the  integral.  i.  fellow;: 


r-l 

)  >/*£ 

»  i 

'  VW 

>rV  >rAr  +  <*" 

i 

^•J 

rr.  -0 

i^-1 

(0*  ’ 

L 

)„  K 

r.4  ^  r*J 

’  il*  £ 
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r.=0 

F-l 

-  y  [(&’  ,,(Iij2))p_-(®D1)KAc  *  ^G' ' '  ^r,2^p-m-l^Dl\r;+lBm] 


sin  0 


:3n(lT-/C)a+\)(XfK 

«/I 7  V l+lt/3  75c- 

*uC(l+t  3)1/2(i+t) 


rHuai+c  2)1/*(i+o'] 

“bi  -  ■  I.  VI  c5  j 
*  >  -  * 


A1  -Vi  (HC/3)1/2c|J 


,-f 'C  2))(l+:«C) 

1  '  -J  7  7 

►  / \  i  •-  -'r,  > 


>- ViZ(l+>C/2)"  ^  .'dC 


P-  +  r^COC.  2)x  IlAiHl 

01  "  -  V7n5  J 


by  re-plating  the  Ej-‘  .  3.41  with  the  eur.r.aticni  ot’ 

Ei;  -  t.  a-.i  its  ••*  :rir.sf  c  *t  v.e  *  err.J  (F1  ’ '(hd)  uni  (G,,,(T|.  ))  , 

if  A  J 

ctuhn  *hf  ft  1 1  :vi r.  g  *_  iua‘ iv  :.n  ft  r  the  rase  n  -  i  , 


(bj(0)r  „ 

^  *  (S->, 


(7(t))r  „ 

IF'  '  (r,  ) )  (K*  )  A  4  (a*"(r,,))  (R"  ).  A  - - =—  -  (3.  ) 

u  r  c  .■  -  f  rti  -  0  VZtf  A*  r 

where-  'c  (t)  ur.i  a.  (*.)  are  giver.  'ey  siting  n  =  1  in  ‘he  first  *w:  cl’  E  jc .  J.J, 
1  > 

f:r  *C  d  <  .  by  E.. ,  f:r  tC.  u  ^ 
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or 


The  surra 


vt)Bl^p 


~,(0  =  \  [si:.  *.  -  4  sir.  Jej 
--  £  [sit.  *,  -t  i  sit. 
bx(0  --•  0,  3j(r)  0. 

ictii,  (bg,)^  at 

y  <VW 


F-i 


S.  )  ,  in  E-1C  . 

4.JJ  r.a; 

rtx  P 

•s.*?  . 

"(V.)) 

(hl.)  a  (O' 

•4  rs  r.+  a 

[(?■■•(-, j))A)r  * 


F- 

> 

r. 


<f' ••('„))  SJ  >.  •  .,)>.«£)  *  , 

*  F  •-  pj  p-* 


(B  +A  ,  ) 
it.  r»+l 


exp  res. 
ar.d  J. 


0:. ?«  tr.e  value-  cf  dlspla server,  t  potentials  have  been  computed,  the 
iicn  fcr  total  hoop  stress  is  obtained  by  ocnvbiniiut  the  lust  of  Evpa.  5. 
with  the  last  of  Eps.  ~.5I. 


(dIU))F  +  t/ilC  I  (F,,’(v))p(P^  )  A  +  (G"'(r,  J){(h'  ).  A 
F-2 

r 

*  ^  [2-  f  ’  ('•  -■.^D4V+i)^far.*A- 

rr.~-C 

t  N/Zr(<’5’,,,,'l  )).<bL>  ■  ry"(\.))  )rl  3.  . 

^  ~  f  ‘  -  I  J.  “i 


wr»*_  r-. 


and  by  E is .  5. 14  fcr  tC^/u  >  2. 

V*>  -!"*=■  5t>  •  =■ 

The  results  of  the  solution  of  the  at;-, 
and  the  corresponding  equations  for  the  ease  n  - 


M‘)  =  :>  . 

•  :  [•  ;eti  f :  r  t  he  hs  s  {  .  *  r 

are  sn.wr.  in  Fig.  -.4  IT r 


43 


v  -  0  and  y  =  1/3-  The  mode,  n  =  0,  dees  not  contribute  to  the  hoop  stress 
at  the  boundary  in  an  incident  shear  wave. 

4  3  Equations  for  Stress  in  the  Medium 

Consider  the  stress  field  at  some  point  away  from  the  boundary  but 
behind  the  incident  wave  front.  The  stress  at  this  point  consists  only  of  the 
incident  stress  if  the  reflected  wave  has  not  reached  the  point,  and  the 
combination  of  the  reflected  and  incident  stresses  if  the  reflected  wave  has 
reached  the  point.  In  Section  4.1  the  equations  found  in  Section  3.1  and  3.3 
for  an  incident  wave  of  dilatation  were  applied  at  the  boundary  in  order  to 
find  the  functions  F^1*2^)  and  Gn+t'(ri2)-  Each  of  these  is  a  function  of  a 
single  argument  given  by  Eqs .  4.3.  These  arguments  are  in  turn  functions  of  r 
and  t.  In  order  to  find  Fr‘+2(r,  )  and  Gn+2(ti2),  a  particular  value  of  r,  r  =  a, 
was  taken  where  the  stresses  were  known.  Once  these  functions  have  been 
determined,  any  value  of  r  may  be  considered  in  the  arguments,  and  the 
variation  with  t  can  be  directly  computed. 

For  convenience  the  case  n  -  1  will  again  be  considered  to  illustrate 
the  method  of  eemputat ion .  The  general  expressions  for  stress  due  to  the 
diverging  wave  given  by  Eqs.  t.6  through  t.8  may  be  used  to  compute  the  stress 
at  any  point  in  the  medium  with  the  upper  limit  given  by  Eq.  4.5.  This  upper 
lirtit  i*  ii.c;  r.venient ,  however,  because  it  leads  to  values  of  ar.d  t,  less 
than  -a/C,  and  -a/C-.,  respectively ,  ir.  ‘he  range  of  integration.  These  values 
represent  time  for  which  the  reflected  waves,  have  not  reached  the  radius  in 
\j-. V  is  more  ..or.  .anient  to  nse  the  upper  limit  which  corresponds  to 
the  first  value  of  F1 '  ’ ( v  )  and  G" 1  ’  (  v)  which  is  r.o‘  zero.  These  limits  may 
be  obtained  from  the  arguments  t(  and  ^  .  Consider  first  the  integrals  involving 
F  '  '  *  i ^ )  when  a  r.ew  variable  t^  is  introduced 
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Then  t,  iii  -r.e  time  elapse  1  after  the  reflected  wave  reached  the  radius  r, 
1 

ar.d  t  is  the  total  tiir.e  elapsed.  Solving  for  t 


r  a 


Substitute  this  into  the  first  of  Eps.  4.3, 


1 


4.56 


The  upper  liir.it  on  £,  £^,  r.ust  be  chosen  so  that  r,^  becomes  -a/c., .  Then 


t 


1 


h  C1 


4.37 


The  variable  of  integration  £ ,  then,  has  the  upper  and  lower  limits 

Vl 

0  <  £  <  51  =  • 

Ir.  the  same  way,  the  limits  on  the  integrals  involving  G1 1 1  (»),)  ray  be  found 
from  n  The  time  variable  +  may  be  defined  as 

Z  <- 


t 


-  t 


r-a 
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or 


Substituting  this  Into  »he  second  of  Ejs.  -*■?  and  recalling  tnat  £  becomes  u£ 
ir.  these  integrals,  cr.e  cbfain: 


1 


I 


The  upper  limit  cn  £  be  corner; 


*2  = 


ctr 


t2Cl 


and  £  varies  in  the  interval 


^5 


*2°! 

o  <  5  <  c,  =  -f-1  • 


In  the  machine  computations,  where  discrete  tiir.e  steps  are  used, 
the  total  tiir.e  is  given  as  before  by  t  =•  pAt .  Also  the  radius  may  be  written 


r  =  a  +  rAt  C^. 

From  Eps.  4.3>  and  A. $8 

~1  -  (p-n)At 
t"2  =  (p-n/a)At  . 

It.  general  the  variable  £  may  be  written 


where  t  varies  fr:m  zero  to  t.  or  t.  ar.d  may  be  expressed  in  general  as  t  =  mAt. 

1  £ 

.’sing  the  above  substitutions,  Eps.  *.36  and  h.39  become 
rlI  -  (p-n)At  -  a/C^  -  mAt 
t,._  -  (p-n/a)At  -  a/C.,  -  irAt 


0  <  rr.  <  (p-r.) 

fcr  v  and 
0  <  r.  <  (r-n/u) 
in  the  j  «a*  lor.  for  . 

Frcr,  Ejs.  a. 57  ar.i  s.*G  the  upper  limits  cn  tne  F"'(r(  )  and  G*  "(ri-,)  integrals 

’t.jV  t-j  w r  .  * 


U-.,)At  C, 


a+;,  At  C 


£  ,  £  = 


(p-n,u)Ai  Cj 
oi  +  **lA*  C  . 


■Since  f  v 3 r ' e s  witr.  radius,  the  p lan*  i* ies  rT^  ,  R?,  .  R^,  .  ,  R^  ,  -:.i 

plf,  must  te  reccr.p  uteu  a*  each  radius  using  Eps.  -*,10  ar.i  a.  11;,  Once  the  new 

JJu 
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values  of  B  have  been  ccr.pu.ted  for  a  particular  radius,  the  equations  for 
stress  in  the  summation  fore  given  by  Eqs.  4.1$  through  4.1$  cay  be  used  with 


the  appropriate  limits. 

If  the  above  computational  procedure  is  adopted  and  if  the  same 
steps  in  time.  At,  that  were  used  at  the  boundary  are  chosen,  the  set  of  values 
of  the  argument  of  F"  "(r^)  used  in  the  confutations  is  the  s«r.e  set  as  were 
used  when  the  original  computations  were  carried  out  at  the  boundary.  Therefore, 
the  limits  of  integration  fall  on  known  values  of  F1 1 1 (q^) .  This  is  not  the 
case  with  G‘ "(ljg)*  however,  since  the  upper  limit  on  ra,  as  shown  above,  is 
p-n/a.  Then  if  m  is  allowed  to  vary  as  m  ---  0,1,2, ...  none  of  the  values  of  t)0 
will  be  the  same  as  at  the  boundary  (see  expression  for  rto  above).  It  is 
simpler  to  perform  the  integrations  to  known  values  of  r\^  (the  same  values  that 
were  used  to  find  the  function  at  the  boundary)  realising  that  these  integrals 
involving  G* ' ' (r^) ,  H2(r,t),  do  not  correspond  to  the  same  times  as  the  integrals 
involving  F‘ ''(q^),  H^r.t).  After  the  integrations  are  performed  and  the 
integrals  are  to  be  combined  in  order  to  obtain  the  stress,  H^(r,t)  must  be 
chi ft el  in  time  so  that  values  of  Hi(r,t)  and  K.(r,t)  corresponding  to  the  same 
total  time  are  combined.  This  shift  will  in  general  involve  interpolation  in 
H_(r,t).  A  parabolic  interpolation  was  use u  for  this  purpose. 

tS 

Ihe  general  expressions  for  .-tress  due  +  the  incident  wave  given  by 
Eqs.  4.2$  uni  4.24  are  valid  at  any  radius  except  that,  u  mtet  be  redefined 
(vi).  Frcr.  Fig.  4.; 


Cjt  -  a  - 
9^  =  cos 


r  cos 
a  -  tC 


( 


r 


1 


-) 


or  for  machine  use 

a,  = 


cos 


,  a  -  pAtC 

( - i) 

'a  +  nAtC, 1 


4.4$ 
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Considering  the  changes  in  and  in  the  integration  limits,  one  has 
the  same  basic  equations  for  stress  in  the  medium  as  were  used  at  the  boundary. 
Tne  total  stress  is  the  sum  of  the  stress  in  the  diverging  wave  given  by 
E-jo  4.15,  4.14,  and  4.15  plus  the  corresponding  stress  ir.  the  incident  wave 
given  by  Ess.  4.25  or  4.24. 


p-n-1 


rr 


j(a,  (t))  +  )  [(P"  1  (ij  )  (Rf,)  A  +  (F’  "(ip.))  .(R^..)  ,B  1 

1'  "p-n  4_,  _  A1  m  m  '  1'  p-m-1'  Al'm+1  mj 


p-n/a-i 


ir.-O 


)  f  (G*  1 1  (r,. ))  (R?, )  A  +  (G1 1  1  (h  ))  .  (R®.  )  _B  !]•  cos  9 

s  |_v  '  *2  p-n.  A1  m  m  '  v  ^"p-m-l  Al'm+l  mJJ 

r.=0 


p-n-1 


■r.  {<h(<»p-r.  *  i  ^  I.'1"  '  ’  Wl>A  *  <F"  1 1'1»,..XU] 

ir.-C 


X  <td)p.,(4>A  *  1g"'I'2)W-i(44a]  } sl" 8 

"•=0 


4.44 


p-n-1 


dd 


•if  i.(t))  +  y  VZc  I  (F‘ "(rtl))  (R„,)A  +  (F' "(h-))  ,(Rnn)  ,B  1 

l  p-n  [_,  |_  v,l"p-r.  Dl'm  m  u/'p-m-l'  Dl'm+1  mj 

r.=0 


p-n/a-l 


Z  ^  f^','(\J)p-r/4)r.Ar.  +  ^"’(^W-I^DI^a]  }  cos  0 
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V.  ESRIVAIIOrf  OP  THE  SHORT  TIME  SOLUTIOfl 

5.1  Eh or t -Tice  Solution  for  Keep*  Stress  at  the  To unit  ry 

Tae  short-time  solution  to  be  described  is  useful  for  several 
reasons:  (l)  it  serves  as  a  check  on  the  machine  solution  of  the  problem  for 
a  short  time  after  tut*  plane  wave  reaches  the  opening;  (2)  it  provides  the 
initial  values  of  the  function  ^+2(  r^)  and  Grt+''(i)2);  and  (3)  it  determines 
the  analytic  character  of  the  solution  at  short  times,  and  gave  insights  useful 
in  developing  the  integration  process.  The  functions  Fn+2(Ti1)  and  Gn+2(r)2) 
which  were  obtained  in  Chapter  IV  by  numerical  integration  are  now  found  by  a 
different  technique,  using  the  same  basic  integral  equations.  Therefore  the 
method  does  not  check  the  derivation  of  these  equations;  it  checks  only  the 
solution  of  the  equations  for  short  time. 

The  same  boundary  equations  given  by  Eqs.  3-41  for  the  case  of  an 
incident  wave  of  dilatation  are  again  used.  The  functions  E^rt+2(T)1)  and 
G  (rj,,)  are  expanded  in  Taylor  series  (with  'unknown  coefficients)  about  the 
initial  values  of  r,  ar.i  t,_  .  Then  the  integrand  of  each  integral  as  well  as 
the  incident  stresses  are  expanded  in  Taylor  series  and  the  integrations 
performed  By  satisfying  the  boundary  equations  for  each  power  of  t  it  is 
possible  to  evaluate  the  unknown  coefficient;,  in  the  series  for  5v*+*"(r|  )  arid 
Gn*“Ci,)  f-uccesil/cly.  Oh'o  the  Series  for  ar.i  Cr‘+,*(r,,)  are  known, 

these  functions  may  be  used  to  find  the  hoop  stn..-;,  yt  tr.v  boundary  or  tli»- 
stresses  in  the  medium. 

Consider  Eqs.  4.6,  4.7  and  4.8  in  which  the  stresses  in  the  diverging 
wave  are  written  in  terms  of  (  in  a  form  suitable  for  combination  with  an 
incident  wave  of  dilatation  for  the  case  n  =  1 . 
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It'  the  value  of  stress  Is  required  at  some  time,  t,  then  the  limits 
_f  integration  of  the  F‘ and  0‘ 1  1  (t|^)  terms  extend  from 

0<£<X  where  £  =  tC^/a  .  5.1 

However,  £  may  be  made  to  vary  from  xero  to  its  upper  limit  by  letting 

t-Pl  .  0<ptl.  5.9 

Then  £  is  a  constant  and  the  variable  of  integration  becomes  p  with  the  limits 
4ero  to  one.  In  order  to  expand  the  function  F1 1 1 (t^)  in  a  Taylor  series, 
consider  the  argument  t|  (Ep.  4. 5), 
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-H1 


o 


a_ 

C, 


£  • 
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wt  r 


°rr  =  '  c3  L  j~b 

i 


a  w  25  TT  145  q2tt2  133  \ 

F  (-  ^-)(5  +  —  PS  +  -52  P  5  +  125  P  5  +  *••) 


4  Aj  |-  (5l(l-p)  -  ^  PC2(1-P)  +  P2I3(1-P)  +  .*•> 

+  A,  tg-)^(3T“(.i-P)3  +  X  P?5  (1-P)2  +  •••)  +  A3(|-)3(3r3(l-3)3  +  ...)] 


dp 


-4^=  f  7k  [a'"(-  §-)(H0  +  +  JaW3  +  ...) 

u^c  ^  Jo  v/2p  *-  C2  5 


2 

B,  (4prd-B)  +  5ap^3(l-P)  +  •••)  +  b2(|-)  (^(l-P) 


if  c£3b(i-p)2  *  ...)  +  b5(|-)';(C3(i-P)3  +  •••)] 


dp 
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,„■&=  r-i. 

*  :3  J35 


f"  1  (-  ^-)(w  +  5p“rJ  +  *  p3^  +  •  •  •) 


*  A,  §-  Upffi-p)  *  5P^J(1-P)  *  ..•)  +  A2(^-)>PC3(1-P)2  +  ...)] 


dp 
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.-.j.r  -.i  n  -i  -  citable  grouping  of  teres 
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also  be  written  in  ’err.:  of  j : wers  :  f  £  by  using  the  series  forr.  of  the 
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trigonometric  functions.  For  n  =  1  the  incident  stresses  are  given  by  Eq.  4.23. 
In  series  fore  these  become 


^(4)  = 


oJK 

3* 


\  25  T  _  2C 

6  ~  S  C  +  3 


5.12 

5-13 

5-14 


In  order  to  satisfy  the  conditions  of  a  traction  free  boundary,  the 
sum  of  corresponding  boundary  tractions  for  the  incident  and  diverging  waves 
must  equal  zero.  The  radial  and  shear  stresses  due  to  the  n  =  1  terms  of  the 
diverging  wave  are  given  by  Eqs.  5-9  and  5*10,  respectively.  These  stresses 
are  given  by  Eqs.  5.12  and  5-13  for  the  incident  wave.  A  series  for  the  total 
radial  stress  at  the  boundary  results  from  the  sum  of  Eqs.  5.9  and  5. 12.  A 
similar  series  for  shear  stress  at  the  boundary  results  from  the  sum  of  Eqs.  5. 10 
and  5.13.  If  these  boundary  stresses  are  to  equal  zero,  each  term  of  the  series 
must  equal  zero.  In  this  way  the  unknown  constants  in  Eqs.  5-9  and  5-10  may  be 


By  substituting  trese  quantities  into  Eqs.  . . }  ar.i  5.4,  cne  obtains 
expressions  for  and  G ‘ 
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a't,(T,5) 


—  a  |  -0.2122  +  0.4421  £  -  0-595^'  +  0,17l6£5  + 

-j7  a  [-0.21WI  +  0.259-12  -  0 . 59251''  +  ...] 


Finally,  x  power  series  for  the  heap  stress  in  the  diverging  wave  is 
found  by  substituting  the  above  constants  into  So .  5.11.  The  total  stress  at 
any  point  is  the  sum  of  the  incident  and  diverging  stresses,  that  is  the  suit  of 
Eps.  5-11  ar.i  5.14. 

°aO  =  0  ^[-1.06701  +  0.9574I2  -  0.2754?5  +  ...  1  5-16 


In  Fig.  5-1  the  results  of  the  short  time  solution  given  by  each  of  Eps.  %  It’ 
through  5.20  for  the  first  five  modes  are  plotted  along  with  the  first  three 
nodes  as  obtained  in  the  computer  solution.  In  the  short  time  solut ion  all 
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—  —7/2 

terns  in  5  out  to  £ were  carried  in  the  computations.  From  these  curves  it 
may  be  seen  that  a  larger  number  of  terms  of  the  series  must  be  used  for  higher 
tnau  for  lower  values  of  n  to  obtain  comparable  accuracy. 

From  the  shape  of  F’ 1 1 (i^)  as  shown  in  Fig.  4.1  agreement  would  not 
be  expected  beyond  one  time  of  transit  in  view  of  the  singularity  at  X  =  2. 

On  Fig.  5.1  the  sum  of  the  first  three  modes  for  agg  as  computed  in  the  machine 
is  compared  with  the  sum  of  the  first  three  and  the  first  five  terms  as 
computed  by  the  short  time  method.  Also  shown  on  Fig.  5.1  is  the  first  three 
modes  of  the  outgoing  wave  combined  with  the  exact  representation  for  the 
incident  stresses.  These  stresses  correspond  to  an  angle  around  the  opening 
of  tJ  =  0°,  cos  nd  =  1.0. 

The  procedure  just  described  is  quite  similar  in  the  case  of  an 
incident  shear  wave.  The  mode  n  =  0  does  not  produce  a  hoop  stress.  This  may 
be  ohown  by  observing  that  the  first  of  Eqs.  3-13  for  bn(t)  is  always  zero  for 
n  -  0.  By  placing  this  result  into  the  first  of  the  boundary  equations 
(Eq.  3 .42),  one  observes  that  F ' ' ( 4 ^ )  is  always  zero  for  this  case.  Then  the 
hoop  stress  in  the  outgoing  wave  given  by  the  last  of  Eqs.  3>39  is  always  zero 
since  F"h,)  -  0  in  the  first  integral  and  sinh  nu  -  0  in  the  second.  The 
dimension less  time,  £,  used  in  the  incident  shear  computations  is  equal  to 
* C .  'a  rutjv.r  char.  tC,'a  as  used  in  the  case  of  an  incident  wave  of  dilatation. 
Tr.e  series  forr.  of  the  hoop  stress  in  the  diverging  wave  only,  has  teen  computed 
fir  r,  =  1  ar.d  These  stresses  are, 


f-  r  . .  -  1 


Ggj  =  a  ^-0. 8294C  +  4.5799?‘ 


1.0197C5  +...] 


5 -01 


for  r. 


ae*  =  ° 


k  -5-0563 

2  * 


8£  +  13.632-r  - 


•75' 
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The  results  of  Eqs .  5-21  and  5-22  are  shown  in  Pig.  5.2  with  the 
corresponding  results  of  the  computer  solution. 

5.3  Short  Time  Solution  for  Stresses  in  the  Medium 

In  Section  4.3  the  method  of  computation  used  in  the  computer 
solution  for  stresses  away  from  the  boundary  was  discussed  for  the  case  n  =  1. 
The  short  time  solution  derived  in  Section  5.1  will  now  be  extended  to  obtain 
the  stresses  in  the  medium.  As  an  example  of  the  computation  procedure,  the 
hoop  stress  away  from  the  boundary  in  the  diverging  wave  will  be  considered. 
The  equation  for  this  stress  is  given  in  general  form  by  Eq.  3*55-  Ety  setting 
n  =  1,  substituting  1  +  C  for  coshu  as  given  by  Eq.  4.2,  and  using  the  limits 
defined  by  Eqs.  4.1,  one  obtains  the  general  expression  for  hoop  stress 


It  is  convenient  in  this  case  not  to  use  the  dimensionless  time  variable  £ 
until  the  last  step  of  the  derivation. 

It  is  apparent  that  Eqs.  4.3  may  be  written  in  the  equivalent  form 

\-1-  (1  *  -r>  “d  (1  * -t4  5-21* 

where  the  dummy  t  ranges  from  zero  to  t.  First  consider  the  argument  .  This 
argument  may  be  written  in  the  equivalent  form 


where 
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The  tine  t^  represents  the  tine  after  the  reflected  wave  has  reached  the  radius 
r,  and  t  varies  from  zero  to  t^ .  This  variation  in  t  may  he  accomplished  hy 
setting 

t  =  pt x  ,  o<p<i.  5.25 

By  using  this  substitution,  may  be  written 
tlx  =  t-^l-p)  -  a/ci  . 

Prom  this  form  of  it  is  seen  that  F' ' 1 (q^)  may  be  expanded  about  -a/c^  in 
powers  of  t-^l-p) .  Therefore,  by  using  the  coefficients  given  in  the  first  of 
Eqs.  5.15,  P"'(tj1)  may  be  written 

(y>  /-  q 

F' ' '  (il1)  =  —  a  -  0.2122  +  0.4421  ~  t^l-p) 

c2  c3 

-0.5956  t2  (1-p)2  +  0.1716  -A  t3  (1-p)3  +  . . .  I  5.26 

a.  fl  J 

In  the  remaining  factors  of  the  integrals  involving  F' 1 1 (q^)  the  variable  £  may 
be  written  as  follows  using  the  substitution  of  Eq.  5.25 

,  TC1  &1C1 

*  ”  r  r 

In  view  of  this  substitution  for  £  in  the  fifth  series  of  Eqs.  5*5  and  of 

F 

Eq.  s.20,  the  first  integral  of  Eq.  5-25,  °qq>  nay  "be  written  in  series  form: 

C„ 


A  C3 


JQQ 


a  f  [-0,2122  +  °'^21  t1(i-p) 

q3  — 

-  0.5956  -A  t2  (i-p)2  +  0.1716  -A  t3  (1-p)3  +  ... 

.  P^C  3  "I  I  t  G 


5.27 


dp 
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This  solution  is  restricted  to  the  case  y  =  l/3  since  this  value  was  used  in 
the  series  expansions  of  Eqs .  5.5  and  the  computation  of  the  coefficients  used 
in  Eq.  5*2b. 


In  a  similar  manner  the  series  expansion  for  G’ 1 1  (t^)  may  he 
determined.  From  Eq.  5-24  the  argument  r\0  may  he  written 


where 


r-a 


The  time  tg  represents  the  time  after  the  reflected  wave  has  reached  the 
radius  r,  and  t  varies  from  zero  to  t0.  This  variation  may  he  accomplished 
hy  setting 

T  =  P  t2,  0  <  P  <  1  .  5.29 

The  argument  r\2  becomes 

q2  =  t2(l-p)  -  a/Cg 

and  the  function  G’ ' 1 (t]2)  may  be  expanded  about  -a/C2  in  powers  of  tg(l-p)  by 
using  the  coefficients  given  by  the  second  of  Eqs.  5-15- 


oc  r  c  c 

o"'(»i2)  =  ~  I -0.2150 t2(i-p)  +  0.2594  -|t®  (l-p)2 

-  0.3989  —  (i-P)3  +  ...J  . 


5.30 


The  variable  £  in  the  integrals  involving  O’  ’ '  (q2)  of  Eq.  5.23  may  be  written 
in  the  following  form, 
tC„ 


By  using  Eq.  5-29  to  vary  t  from  zero  to  tg  and,  by  letting  =  aC^,  one  may 
rewrite  this  expression 
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£  = 


e  t2  ac1 


It  follows  from  Eq.  5-50  and  the  above  substitution  for  t,  in  the  sixth  series 

Q 

of  Eqs .  5-5  that  the  second  integral,  o0e  ,  of  Eq.  5.23  may  be  written  in 
series  form. 


1 


G 

aee  -  0 


p  (J  Q” 

J  ~  j  -0.2150  t2(l-p)  +  0.2594  t2  (1-p) 


0.3939  x  %  +  •••!  {-#-  /  1 

2  J  ^  a  V  2aptgC1 


40 


apt2Ci 


apt-C,  2  „  apt  c  3  n  .  t„ac, 

t  d  -M  .  L  / _ ?..A\ 


o(-^)  +§(- 


1)  + 


)  dp 
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The  following  equations  result  from  performing  the  integrations  on  the  right 
sides  of  Eq.  5.27  and  5.31,  for  the  case  \  =  p,  v  =  1/4,  yielding  a  =  -I _ 

yr 


F 

=  a 


t  c  n  v  v 

i-i  +0.4244  -  0.5894  t,  -r  +  0.4199  t2  1 


*1  a 


1  a2 


C3  Q 

0.1569  t l  -I  -  0.4598  tx  —  +  0.3832  1 

a'' 


ar 


c?  „  c?  ,  c? 


0.1949  tj  -i-  -  0.4377  t2  -I  +  0.2605  — 


a  r 


ar 


0.0^97**  t 


3  ^i+  1 

1  „3  •••'] 


5.32 


Jee 


[-0.5227  ^+0.36041! 

-  0.1617  T2  ...  ]  . 


a2r 


5-33 


Tne  results  obtained  in  Eqs.  5 • 32  and  5*33  may  not  be  added  directly.  Equal 
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values  of  and  tg  correspond  to  different  times  t,  being  the  interval 
after  the  reflected  dilatation  wave  has  reached  r  and  tg  being  the  interval 
after  the  shear  wave  has  reached  the  radius  r.  For  this  reason  the  two  parts 
of  the  hoop  stress  will  be  evaluated  and  compared  with  the  machine 
solution  separately.  Equations  5-32  and  5-33  may  be  written  in  terms  of  the 
dimensionless  times 


Vl 


and 


as  follows 

+  0.4244  -  0.5894^  +  0.4199^ 

-  0. 15695 l  -  0.4598?!  7  +  0.3832?®  |  -  0.1949?^  7 

2  2  3  -1 

-  0.4377?!  Sg  +  0.2605 ~2  ~  o.o«*974?J  ~  +  ... 


5.3^ 


'08 


=  «  [-0.5227?!  7  +  0.3604?^  |  -  0.l6l7?|  4  +  •••]* 


5-35 


In  Fig.  5-3  the  stresses  given  by  Eqs.  5-34  and  5-35  are  shown  along 

with  the  corresponding  results  from  the  computer  solution  for  a  value  of 

—  =  1.25.  It  may  be  seen  from  these  curves  that  the  difference  in  the  two 
a 

solutions  is  quite  small  out  to  ?!  and  ?g  =  0.6. 
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VI.  RESULTS  AND  H! TERPRETATION 

o.l  Incident  Wave  of  Dilatation 

The  results  of  the  computer  solution  for  hoop  stress  at  the  ■boundary- 
due  to  three  terms  of  the  Fourier  series  representation  of  an  incident  wave 
of  dilatation  are  shown  in  Figs.  6.1  and  6.2.  The  most  important  information 
which  follows  from  an  examination  of  these  curves  is  summarized  in  the 
following  table. 


V 

Q 

Maximum 

Stress 

°ee/° 

Time  of 

Max.  Stress 
Transit  Times 

Static 

Value 

°Ja 

0 

0° 

+1.29 

4.8 

+1.00 

90° 

-3-33 

3-3 

-3.00 

1/3 

0° 

+0.11 

4.0 

0 

90° 

-2.93 

3-7 

-2.67 

In  Sections  5-1  and  5-2  a  method  was  discussed  which  provides  a 
check  on  the  machine  solution  of  the  basic  equations  at  short  time.  After 
several  transit  times  have  passed  the  stresses  should  approach  the  static 
values  as  given  above.  This  provides  not  only  another  check  on  the  solution 
of  the  basic  equations  but  to  some  extent,  on  the  formulation  of  the  solution 
to  the  problem.  The  static  stresses  in  the  case  of  plane  stress  are  given  by 
the  well  known  Kirsch  formulas,  which  are  shown  as  follows  in  Ref.  (ll): 

arr  =  I  (1  '  +  I  (1  +  '  ^2>  COS  2® 

r  r  r 

Qqq  -  %  (l  +  ~  ^  (1  +  cos  20  6. 1 

r  r 

ars  =  - 1 (1  -  ^ +  ^  8in  20  * 

r  r 
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These  equations  may  be  used  to  solve  the  plane  strain  problem  by 
superimposing  on  the  stress  field  a  the  stresses  due  to  another  stress  field 
acting  in  the  perpendicular  direction  and  of  magnitude  v  o.  By  this  means 
the  static  solution  for  agQ  around  the  opening  as  well  as  stresses  in  the 
medium  can  be  found  for  all  the  problems  considered,  and  in  all  cases  the 
stresses  approach  the  static  values  quite  closely  after  5  to  7  transit  times. 

The  results  shown  in  Figs.  6.1  and  6.2  are  compared  in  Figs.  6.J  and 
6.4  with  the  corresponding  results  from  Refs.  (3)  and  (4).  The  hoop  stress  on 
the  boundary  at  0  =  0°  from  Ref.  (3)  is  shown  in  Fig.  6.3  to  be  quite  different 
from  that  given  by  Ref.  (4)  and  the  present  work.  Also  the  hoop  stress  at 
0  =  90°  from  Ref.  (3)  is  shown  in  Fig.  6.4  to  be  less  than  that  given  by  Ref.  (4), 
but  the  difference  is  not  as  large  as  for  0=0°.  It  is  interesting  to  note 
that  all  three  computations  give  quite  close  results  for  the  maximum  values 
of  tensile  and  compressive  stresses.  In  Ref.  (4)  the  reason  for  the  difference 
in  stresses  given  by  Ref.  (3)  is  explained.  In  using  the  Duhamel  integral  for 
computation  of  the  boundary  tractions  from  the  diverging  step  wave,  singularities 
of  the  same  type  as  in  the  present  work  occur  at  values  of  tC^/a  =  0  and  2.0. 

For  this  reason  great  care  must  be  taken  in  the  numerical  integrations  in  this 
region.  In  Ref.  (3)  a  constant  time  interval  was  used  throughout  this  integra¬ 
tion;  however,  as  shown  in  Ref.  (4)  it  is  necessary  to  reduce  the  length  of  time 
interval  in  the  region  of  the  singularities.  It  is  suggested  that  the  integra¬ 
tion  technique  discussed  in  Appendix  A  would  have  been  useful  for  the  computa¬ 
tions  of  Refs.  (3)  and  (4)  since  the  singularities  have  the  general  form  of 
-l/2 

t  '  .  The  method  discussed  in  Appendix  A  is  a  more  accurate  means  of 

accounting  for  the  singularities  than  is  provided  by  reducing  the  length  of  the 
integration  interval. 

The  hoop  stresses  at  0  =  0°  and  90°  from  Ref.  (4)  agree  with  the 
present  work  as  closely  as  the  curves  in  Ref.  (4)  can  be  read,  except  during 
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the  first  transit  time.  With  regard  to  this  difference  it  should  be  recalled 
that  the  incident  stress  was  represented  by  three  terms  of  Fourier  series  and 
that  diverging  wave  stresses  were  computed  to  correspond  to  these  three  terms. 
Once  the  diverging  wave  stresses  are  found,  they  must  be  combined  with  the 
incident  stresses  in  order  to  find  the  total  stresses.  However,  it  must  be 
decided  whether  the  three  modes  of  diverging  wave  stresses  should  be  combined 
with  the  same  three  modes  of  the  incident  wave  stresses  or  with  the  total 
incident  stresses  which  may  be  represented  exactly.  In  Figs.  6.3  and  6.4  the 
results  of  both  methods  of  representing  the  incident  stresses  are  shown.  Both 
Rex's.  (3)  and  (4)  apparently  used  the  exact  representation  of  the  incident 
hoop  stress;  however,  the  hoop  stresses  shown  in  Figs.  6.1  and  6.2  as  well  as 
the  stresses  in  the  medium  and  the  stresses  caused  by  decaying  waves  found  in 
the  present  investigation  were  computed  using  the  first  three  modes  of  the 
series  representation  for  the  incident  stresses.  It  was  not  possible  to 
determine  the  cause  of  the  disagreement  between  the  results  for  hoop  stress 
1‘rcm  Ref.  (4)  and  the  present  work  during  the  first  transit  time  when  the  same 
incident  stress  was  used.  In  view  of  the  excellent  agreement  with  the  inde¬ 
pendently  computed  short-time  stresses,  the  present  computer  solution  acquires 
considerable  credibility  at  least  for  early  times.  Moreover,  if  one  considers 
the  presen*  method  of  solution,  it  hardly  seems  possible  that  the  hoop  stress 
could  be  considerably  in  trror  during  the  first  quarter  transit  time  without 
causing  a  large  difference  at  one  transit  time. 

Th^  use  of  the  total  incoming  wave  causes  the  total  hoop  stress  at 

0  -  0U  to  have  an  improper  initial  value.  This  may  be  shown  in  the  following 

way.  For  y  -  1/5,  a  -  0°  and  the  incident  stress  a  -  -1.0,  the  incident  hoop 

stress  is  \rs  =  -1 /?;  however,  at  the  same  point  the  diverging  radial  stress  must 

equal  a  =  a  -  +1  at  t  =  O  and  therefore  the  diverging  hoop  stress  aOQ^va  =+ 
rr  yy  rr 
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since  it  may  be  argued  that  there  is  no  strain  in  the  9  direction  for  9  =  0  at 
this  time.  Therefore  the  diverging  wave  should  have  an  initial  value  of  hoop 
stress  equal  to  +1/3  if  all  the  modes  were  considered,  and  the  initial  value  of 
total  incident  and  diverging  stress  at  9=0°  combined  should  be  zero. 

In  Figs.  6.3  and  6.4  the  hoop  stresses  at  9  -  0°  and  $/Q°  are  shown 
using  the  first  5  modes  as  given  in  Ref.  (4).  By  comparison  of  these  curves 
with  those  given  for  3  modes,  the  magnitude  of  error  may  be  estimated.  The 
error  appears  to  be  quite  large  at  9  =  0°  for  the  first  two  transit  times  as 
would  be  expected.  At  90°  the  addition  of  two  more  modes  makes  very  little 
difference.  Also  the  use  of  5  modes  gives  the  initial  value  of  at  0  =  0° 
much  closer  to  the  required  value  of  zero. 

By  comparing  Figs.  6.1  and  6.2  it  may  be  seen  that  a  reduction  in 
Poisson's  ratio  increases  both  the  maximum  tensile  and  maximum  compressive 
stresses. 

In  Figs.  6.5  through  6.8  is  shown  the  variation  of  stress  with  time 
at  several  radii  for  v  -  0  and  for  v  =  1/3  due  to  an  incident  wave  of  dilatation. 
These  curves  show  the  free  field  stress  ur.tU  the  reflected  wave  reaches  the 
point,  at  which  time  there  is  a  sharp  change.  At  some  time  after  the  reflected 
wave  passes  the  point,  however,  the  stress  approaches  the  long  time  or  static 
value  consistent  with  the  particular  point.  From  these  stress-time  plots  at 
various  radii  it  may  be  seen  that  the  dynamic  effect  on  the  peak  stress  is  not 
as  great  away  from  the  boundary  as  at  the  edge  of  the  opening.  The  static 
variation  of  stress  with  radius  as  found  from  Eqs.  6.  1  are  shown  in  Fig.  6.  9. 
Also  shown  in  Fig.  6.  9  is  the  variation  of  radial  and  hoop  stresses  with  radius 
at  several  times. 

In  Fig.  6.7,  which  shows  the  radial  stress  Vs.  time  curve  for  9=0° 
and  v  =  1/3,  an  indication  of  the  mechanism  which  might  cause  spalling  is 
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observed.  The  stresses  should  be  considered  only  qualitatively  during  the 
transit  of  the  wave  across  the  opening  since  only  three  modes  are  used.  At 
r/a  =  1.25  for  example  the  stress  drops  very  sharply  from  -1.00  to  -0.05  when 
the  reflected  tension  reaches  the  point.  If  the  incident  stress  wave  had  a 
decay  behind  the  front,  this  drop  in  stress  could  easily  produce  a  total 
tensile  stress  A  sufficiently  large  tension  would  cause  failure  of  the 

medium.  Since  the  slab  bounded  by  the  region  of  tensile  failure  and  the  free 
surface  of  the  opening  has  a  particle  velocity  equal  to  the  sum  of  that  in 
the  incident  compression  wave  plus  that  in  the  reflected  tension  wave  (both 
in  the  inward  direction)  it  would  tend  to  fly  inward.  This  slab  remains  con¬ 
tinuous  with  the  adjacent  material  in  a  circumferential  direction,  however, 
and  the  energy  trapped  in  it  must  be  sufficient  to  break  it  away  to  form  a 
spall. 

From  the  results  of  an  incident  step  wave  it  is  possible  to  find 
tne  stresses  caused  by  an  incident  wave  of  any  other  form  using  the  Duhamel 
integral. 

dmn  =  %nIS)(l)-/p  ^  "  T>]  dT  6.2 

where  av  is  any  stress  due  to  a  unit  step  wave,  and  a  is  the  corresponding 
niii  inn 

stress  du»  to  an  arbitrary  stress  wave  F(.t)  To  illustrate  the  use  of  the 
Duhamel  integral  and  study  the  effect  of  a  decaying  stress  wave  the  hoop  stress 
was  computed  at  «  =.  O'  and  90°  for  v  =  l/3  and  a  =  17-5  ft.  The  variation  of 
stress  behind  the  wave  Trent  was  assumed  to  be  of  the  following  form 

P(t)  =  Pc(l  -  f-)  e“Kt/t° 

~o 

where  K  is  a  constant  and  is  the  total  positive  duration  of  the  wave.  Hoop 

stresses  are  shewn  in  Fig.  6.10  and  6.11  for  P  =  4CG0  psi,  t  a  1.8  sec., 

o  o  1 
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Cx  =  17,300  ft. /sec?,  and  K  -  ICO,  2J0  and  4C0.  Also  shown  is  the  stress  due 
to  a  step  wave  with  a  stress  of  40C0  psi.  By  increasing  K,  the  rate  of  decay 
behind  the  shock  front  is  increased.  The  values  of  K  used  represent  a  rather 
sharp  decay  which  wculd  be  consistent  with  high  stresses  produced  by  a  large- 
yield  nuclear  weapon.  The  peak  compress ive  stress  at  0  =  §0°  is  reached 
within  5  transit  times,  or,  in  terms  of  the  wave  velocity,  only  the  first  6  ms 
of  the  wave.  In  the  case  of  the  tensile  stress  the  peak  occurs  within  5  transit 
times  or  within  10  ms.  after  the  wave  reaches  the  opening.  Therefore  only  the 
very  early  portion  of  the  wave  history  is  pertinent  in  determining  the  maximum 
stress.  In  5  ms.  the  pressure  lias  decreased  to  5020  psi  for  K  =  100  and  to 
1516  psi  for  K  =  400. 

In  general  only  the  maximum  tensile  and  compressive  stresses  have 
been  shown  since  these  stresses  are  critical  in  determining  the  fhilure  of  the 
medium.  The  computer  code  which  was  written  to  compute  these  stresses  can  also 
compute  stresses  at  any  number  of  evenly  spaced  angles  around  the  opening  and 
at  any  radius  out  from  the  boundary  so  long  as  the  radius  is  an  integral  number 
of  distances  AtC^.  This  code  is  considerably  less  complex  than  that  used  by 
the  authors  of  Refs.  (5)  and  (4).  To  solve  a  problem  in  which  the  stresses 
are  computed  at  15  equiangular  increments  around  half  the  opening  and  for 
sufficient  time  to  reach  the  static  solution  requires  approximately  20  minutes 
of  machine  time  for  5  modes.  For  a  problem  in  which  the  stresses  are  computed 
at  l/t  radius  spacing  for  4  radii  cut  from  the  bcar.aary  and  with  the  same 
angular  increments  requires  approximately  >0  minutes. 

It  is  of  particular  interest  to  notice  the  effect  of  changing  the 
length  of  time  steps  in  the  computations.  The  following  table  shows  the  hoop 
stress  at  the  boundary  for  0  =  90°  and  v  =  l/j  at  several  times  as  the  length 
of  time  steps  are  varied. 
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lire  _,*ep.;, 
per  transit 
t  ime 

I  ri(:  »reos 
t>  ,/cs 

*W «  at 

tC^/a-1.0 

aed/a  at 
tC^/a-2.0 

.-0 

-0.7^56 

-I.83&7 

-  - .  yo  - 

-U.79U-, 

-1.8334 

■j0 

-0 . 7652 

-i.8336 

It  i.i 

-4.d*0h 

-O.7P29 

-I.8336 

The  effect  ef  ^ 

the  . e.igtn  _f  j 

Live  step  on  the  peak  stress  is  negligible. 

At  earlier  times  the 

fleet  is  larger, 

but  still  unimportant  from  a  practical 

standpoint.  Tne  value 

of  40  steps  per 

transit  tine  was 

used  for  the  computation 

of  stresses  iu*.  to  an 

incident  wave  of 

dilatation.  The 

incident  shear  problem 

was  found  to  be  at  core  ..ens  it ive  to  changes  in  the  length  of  time  step, 

so  £0  steps  per  transit  time  was  u-ei  in  those  computations. 

t.2  Incident  Shear  Wc. e 

In  figs,  o.ij  and  6  it  are  shown  the  results  of  the  computer  solution 
for  h.  op  ot  re.  .  cn  ‘he  t '  lury  due  to  three  cede.;  of  an  incident  shear  wave 
vi  h  msgi.ituie  of  str-n...  i.  Tne  ir.p.ortant  information  shown  by  these 

curves  is  s'jr.tarlo*. i  in  the  to) lowing  table. 


:s,.imc.  lire  of 

-th'-  !*u\.  atr..;s, 

c,t  i  i  Transit  Tic13;- 


<  4.  CO 
-4.00 

<•4  .CO 
-4.  CO 


The  rt re •**. *ji. i  *i.e  ■  fti.'ng  due  t  i  shear  wave  ..houH  approach  the  static 
values  alter  a  i>.  .g  tire,  iiies*-  e’uti’  values  r.uy  he  ottasi.ei  using  Eqs.  6.1 
by  s  jj  s*  is.p  ;:.ii  g  ■.  j  .-4 1  ..tpr.;r..v  i.T  ui.j  *e.io  ii .‘less  fields  acting  $QV  apart. 

In  this,  way  i  field  ■  t  |  irt  hear  is  create  1  er..  ied  t he  upening.  The  maximum 


i  /  r 


, .  S'. 


-  i  'ih 


Static 

value 
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static  stresses  obtained  in  this  way  are  found  tc  occur  on  the  boundary  at 
45°  and  135°  measured  from  the  direction  of  wave  travel.  These  static  values 
of  maximum  hcop  stress  are  shown  in  the  above  table. 

Also  shown  in  Fig.  b.ij  is  the  hoop  stress  at  8  =  45°  due  to  an 
incident  shear  wave  as  shown  in  Ref.  (12)  for  v  ^  1/3,  In  Ref.  (12)  the 
authors  of  Ref.  (5)  use  the  same  method  that  they  used  to  solve  the  incident 
dilatational  wave  problem  to  solve  the  incident  shear  problem.  Therefore,  a 
constant  time  step  in  the  Duhamel  integral  near  the  singularities  at  tCg/a  =  0 
and  2  was  probably  used  as  before.  This  would  explain  why  the  curve  of 
Ref.  (3)  does  not  agree  with  that  from  the  present  computation.  The  case  of 
an  incident  shear  wave  was  not  solved  in  Ref.  (4). 

Also  Ref.  (12)  shows  a  hoop  stress  at  6  =  0°  which  is  not  found  by 
the  present  theory  and  which,  in  fact,  can  be  shown  to  be  nonexistent  by  a 
simple  symmetry  argument.  For,  by  inspection 


Therefore 


5u  Jv  , 

7  -  ~  +  —  and  v  are  even  in  y , 

xy  oy  ox 


dv  .  du  , 

and  are  even  in  y  , 


This  requires  that 


du 

u,  =  e  , 
’  ox  x 


3v 

37 


be  odd  in  y  . 


Since  •:  and  f  are  both  odd  in  y,  they  ire  ,*ero  at  y  -  0,  and  as  a  consequence 

X  j 

the  stresses  a ar.d  uyt(  are  sero  at  y  0. 

The  passage  of  a  shear  wave  across  an  opening  results  in  a  higher 


concentration  factor  for  both  tension  and  compression  than  the  passage  of  a 
wave  of  dilatation  and  therefore  may  cause  the  more  critical  situation. 
However,  the  intensity  of  stress  in  a  shear  wave  is  much  less  than  the  stress 
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jl-.  u  vi'-e  cl'  diiutati;:,  l\  r  at.  air  induced  grcur.d  shock,  and  may  be  absent 
altogether  fcr  the  iirect  ground  ^noex.  Tne  effect  on  the  stress  concentration 
factors  of  varying  Poise :a*s  ratio  is  snail  as  shewn  by  comparing  the  curves 
of  fig.  *0  12  Und  o. it. 

The  computer  cede  written  to  solve  the  incident  shear  wave  problem 
for  Keep  stress  at  the  tcardury  is  puite  similar  to  that  discussed  for  an 
incident  dilatational  wave,  ard  therefore  requires  approximately  the  same 
amount  of  time  for  solution  on  the  machine.  This  code  was  found  to  be  some¬ 
what  more  sensitive  to  changes  in  the  time  interval  than  was  that  for  an 
incident  dilatational  wave.  The  reason  for  this  is  probably  the  sharper  slope 

of  the  Fourier  coefficients,  a  t*)  and  b  It;,  in  the  series  for  incident  shear 

ri  n 

stress  when  the  wave  fr:nt  is  near  tne  front  and  reur  of  the  opening. 

t  J  '  -e  ir.  Protective  Construct icr. 

Tne  results  of  thlw  investigation  may  be  used  as  a  guide  in  the 
st.iy  of  effects  . f  „tr*.-..«  wave.,  ot  openings  in  rock.  Among  the  first  things 
to  te  decided  ir.  such  u  study  is  the  type  of  failure  to  be  expected.  For 
example  the  maxlmor.  c  empress  ive  s’ res  s,  which  occurs,  rat  the  ends  of  the 

>  p  r  i.-  1  •  -wa.-e  front  tor  it. client  waves  of  dilatation,  may  cause 

cr-.u-r.,.-  -f  iinr.  2.  *he-e  —  cr.c  Tr.e  ~axitr  ot  tensile  stress' ,  which 

.■csjs.  a-.  i:  ci  Vv  i.  .*.ter  t  erpu.  ii  c  i  far  *  -  the  wave  front,  may  cause 

,-r.  t  V,  l.  it,  r-.—T-.s.  'reTirg  ic-oS'e  reck  which  would  fall  into 

the  opM.is  *.  sir.-f  » •  lie  s*rer.g*h  ef  material.,  sack  as  rock  is  much 

vr,  i;.  '  r  :■  se.  :• .  .•<;  - 1  re  g*h,  termich  failure  may  occur  at  a  free  field 
street  which  i.-  not  sufficient  to  cause  compression  failure.  Thus  the  ratio 
of  maximum  ’ er, s i  1  >•  tv  rssxis.am.  csr.press.ive  stress  becomes  important  in 

deciding  which  thes-  *ype^  of  failure  w.i  c-cur.  Another  type  of  damage 
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and  possible  failure  nay  c  . .  ,r  fr._n  spall i  .g  which  was  described  earlier.  If 
tiie  free  field  stress  is  nigh  enough,  the  materia L  ir  the  entire  region  of  the 
opening  may  fail  ar.i  close  the  opening.  V.iis  ceu.'d  e .•  .nr  if  the  opening  were 
sufficiently  close  to  the  center  of  burst  of  a  nuclear  weapon.  Construction 
of  a  lining  within  the  cavity  would  afford  some  protection  against  spalling 
or  loose  rock. 

Though  the  maximum  stresses  in  the  region  of  a  cylindrical  opening 
have  been  determined  for  a  homogeneous  medium  and  ideal  wave  forms,  there  is 
still  considerable  uncertainty  in  the  practical  case.  The  results  obtained 
may  only  be  used  as  a  guide.  For  example  rock  contains  faults  and/or  joints 
and  therefore  is  often  weaker  than  the  conventional  strength  tests  would 
indicate  Furthermore  rock  may  be  so  cracked  or  so  damaged  by  the  mining 
operation  that  it  has  no  tensile  strength  near  the  opening.  Also  the  wave 
form  is  surely  altered  in  traveling  through  a  nonhomogeneous  medium  and  probably 
dees  net  have  a  sharp  front.  The  effect  of  a  rise  time  would  be  to  reduce  the 
dynamic  stress  concentration  factor  and  to  reduce  the  possibility  of  spalling 
failure  From  this  discussion  it  ray  be  seen  that  considerable  judgment  is 
required  in  making  allowances  for  these  uncertainties. 

In  Fef.  (iji)  many  of  the  practical  problems  encounter’  in  the 
application  of  theoretical  results  of  the  type  obtained  in  the  present  work 
are  discussed  Also  the  various  wave  types  are  discussed  and  a  method  for 
ccmbinirg  *he  effects  of  waves  ■  }’  JiJ -tt  at  ion  n.  i  wav* us  wr.er.  I .  th  are 

created  by  an  air  blast  wave  moving  on  the  surface  is  presented. 
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r<-.--se.  it.  »r.e  viciiity  . f  x  _*y 1  ii.drieui  cavity  due  to  two  types 
v,f  stress  w vec  i  x'.v  own  stuliei.  Aj.v  the  effect  of  changing  the  shape  of 
the  stress  wave  behind  the  front  ana  <_  f  varying  the  properties  of  the  material 
have  tees.  cons  iderei  t  r-..t  the-e  -tui.es  certain  cone  fusions  may  he  made 
concerning  the  r.etnc i  of  solution  use!  to  solve  the  problem  under 


consideration. 

(i)  Tr.e  method  of  solution  iescr.tei  it.  the  present  work  can  be 
usoi  . .ccecsfuily  uni  efficiently  to  solve  the  problem  under  consideration. 
Cr.co  the  -.lotion  for  u  unit  s’ep  p  .  .  o  ■  •  pj  /v. own,  it  is  possible  to  find  the 
^tr*:. ...  for  x  •»  j  ■  q-  :  i.  *  t.  •-  -1  integral . 

(;•}  Tr.-c  —tort  time  .•.elution  may  ce  used  to  check  the  initial 
portion  of  in*,  nuoh.i.e  - :  :ut fy  ... ir.c  u  larger  number  of  terms  in  the 
luyior  e:*.jjt..oi::.o,  tt.u-  ..uj’k:.  could  ce  -..iteniei  to  apply  to  longer  times. 


appear-  to  give  maximum 


v. -.:.u*y  uv.-y  wit-,  .ufficient  accuracy  for  practical 

I:,  rd**r  t.-  . ■*  ujy  *he  .-.pulling  f hencre:. on,  however ,  it  is  necessary 
.  .  ;  :  j  rt  times  are 


.  -.i  or 


.it-  .  •  *  j  - .  .  i  '•  i  .  :  ,r-..-.r  •••iieo,  eltr.er  by 

. .  u  .  . . .  r  ■  j  j.  .v..  .  d  *.  other  problems.  A 

‘  *_*w  *  ■.  i  i  '  ’  !*t. 1  l  •  .i1  ;  -  J  f-  _  x  1  c 

;  **  ■*  t  ,*.jr  ■  r  .  •  ;  i‘  i-. .  'i'll  1  5t.  Appendix  G  could 

vrif*  si.  f ;  r  .«  z.:t\  ur  i  .  ufi'l  i*  ’  *  r  x  tUii:,  ur.y  -accuracy 


a- 
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(2)  By  using  a  large  number  of  xs.co.es  either  with  the  computer 
solution  or  the  short  time  solution  the  stresses  at  early  time  could  be  obtained 
with  sufficient  accuracy  to  make  a  complete  study  of  the  spalling  problem. 

(J>)  A  complete  study  of  various  wave  forms  such  as  effect  of  rise 
time  and  decay  rate  would  be  possible.  If  wave  forms  with  rapid  changes  in 
stress  were  to  be  considered,  it  would  be  necessary  first  to  obtain  the  solution 
with  more  modes  included. 

(4)  The  problem  of  a  cylindrical  opening  with  a  lining  could  be 
solved  by  the  method  presented.  Either  the  case  of  a  thin  lining  or  that  of 
a  thick  lining  could  be  considered. 
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PsfUYAIlC  .  Ui  id2  M-MT-J.^X  FiiifJlA'A",';  >CBTOA 

Fr:-«  the  for-  of  re  ' r  .... .  equations  feu  diverging  wives  written 

in  terra  _n  £,  E*i.  '*.0  -  '«.c-  for  e a- •-.{•!«,  i'.  3.uy  ce  seen  that  one  terr.  of 

l/2 

each  expr-.Mcn  has  or.  it,*  egr  .•*!•.  • '  •  grurl'y  at  £  -  0  due  to  *he  terr  £  '  in 

the  denominator.  It  is  apparent  that  a  trapezoidal  or  parabolic  approximation 
would  not  give  a  finite  result  in  the  f 1  r„.t  _tep  of  the  numerical  integration 
even  though  the  Integral  converges.  Ar.  integration  procedure  was  derived  to 
avoid  this  difficulty,  ur.d  is  describe a  in  this  appendix.  After  this  procedure 
had  been  ir.  use  for  seme  tire,  it  was  found  that  u  very  similar  procedure  was 
described  by  Jeffreys  ar  i  Jeffry,  in  E-*f.  (?)• 

Consider  the  second  integral  involved  in  the  computation  of  a  Q  in 
mode  1  given  by  E p,  A.  7. 


5 

1  iV  '  •  | 

,  .  ril(i,  or,. utf'V 

1,  J 

L-  t  uii  tn*.  integral  1  e.-:  cm  *  1  £  ' "  t"  upproxiruted  by  a  linear 
relation  in  t*  -  Interval  ft  v  <* 


'•’’’(a.) 


[ n( j  <  y-t <buwC,'-t.ir-'r.-’1 

~~F=~ 7====  ;  v  1  ■* 

1  w  i'  f  1+ 0  ore; 

“■  t 


<5-^0 


k  .1 


in  the  interval  5  =  rxif;  to  £  -  (r.+  l)i£.  With  this  substitution,  the  integral 
in  this  interval  becomes, 
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X,, 


'niC  V£  L 


^+l“^ 


(5-cA£)] 


d£ 


[^?Vi  *  (|  tV2)  _  ^  ^(ati/2, j'" 

m££ 


(m+l)&£ 


A. 2 


S’Tjiit’hrforwur-i  evaluation  of  this  expression  leads  to  the  following  value  for 
the  left  side  of  Eq.  A. 2. 


|"g^  ^(ItoJUk+I)1'^-!^2}  -  j  [  (rc+l)^2-m^2]  j- 
♦  6to+1  '(j  C  (m+l)3/fi-m3/'S  )  -  2m[(m+l)l/2-ml/2]}  j 


A. j 


This  expression  gives  the  value  of  the  integral  between  and 
(ir.+  lp£  as  long  as  the  values  of  g^  and  ^,+1  are  defined  according  to  Eq.  A.l. 
Tne  fur  cticns  of  r.  by  which  each  ordinate  is  multiplied  differs  for  each  step, 
since  r.  charges,  but  will  be  the  same  for  integration  of  F1 ’’(q^)  and  G’ 1 1 (q^) 
terns  Also  ,  this  expression  is  not  limited  to  n  -  1;  it  is  a  general 
integration  formula  for  which  g  is  any  function.  From  Eq.  A. 5  let 


A 

r 

B 

c 


-  U+n)  [(r.+  l)1  “-ir.1'-]  -  j  [  (c+1) ] 

y  [(n+l)',/‘i-r^  “  ]  -  ] 

j 


A  k 


HU.  '.nose  at  ore",  iut  ions  the  integral  to  be  evaluated  nay  be  written 


I  • 

J 


Jr.  4  1  PC 


V"(v,) 


•  •  r 


'jUxi  V  laaC/2 


adC 


(G"’(r(,))p.r 


R°  A 

r.  tn 


+  U 


vTl  )) 


c  P 


RU 
•n-1  n+1 


A. 5 


where 


RG  is  the 

r. 


factor 
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u(  l+5at+6acy +20^)3 

•/an  c? 

p 

evaluated  at  £  =  cud;  and  ET+1  is  the  same  factor  evaluated  at  £  -  (rc+1  )£.{; . 

In  computing  the  factors  A  and  B  ,  the  small  difference  between 

C*  Hi 

large  numbers  are  involved  when  m  gets  fairly  large  if  the  computation  is 
made  in  a  straightforward  manner.  Therefore  these  weighting  factors  were  also 
written  as  asymptotic  series  which  were  used  after  m  reached  a  certain  value. 
These  series  are  as  follows 


»  1 

fl  1  .  1  1  ,  1 

A  —  _ 

m  2  •Im 

1  3n  2  3  J 

L  1cm  32r* 

_  1 

",  1  1  1  ] 

B  =  - — 

m  2-i/m 

L  >  16m-  8m'  J 

APPENDIX  B 

SOME  FE0F2RTIES  OF  THE  DISPLACEMENT  POTENTIALS 
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B. 1  General  Solution  of  the  Wave  Equation  in  Cylindrical  Coordinates. 

lr»  Section  3*3  it  is  .shown  that  tvo  hyperbolic  partial  differential 
equations  in  two  independent  variables  result  from  substituting  Fourier  series 
for  <p  ar.d  y  into  the  wave  equations.  These  differential  equations  are  Eqs .  3*21. 
The  relevant  parts  of  the  general  solutions  of  these  equations  were  stated  as 
Eqs.  3.22.  It  will  now  be  shown  that  these  solutions  satisfy  the  differential 
equations.  Since  both  differential  equations  and  both  solutions  are  of  the 
same  form,  it  is  only  necessary  to  show  the  validity  of  one  of  the  solutions. 

Then  it  is  to  be  shown  that  for  n  =  0,  1,  2,  ... 


,  ,1  d,r* 

*n  ”  'r  dr'  c 


B.l 


satisfies  the  equation 


_  d^f  -  ,  df 

P2  / _ n  dn+1  n,  _  -a 

C1  2  r  dr'  ' 

dr 


B.2 


I*-  is  first  assumed  that  f  is  a  solution  to  Eq.  B.2.  Then  it  will  be  shown 

n 

that  f  +  given  by  Eq.  B.l  is  also  a  solution,  and  therefore  the  general 
edition  for  any  value  of  n  is  given  by  Eq.  B.l  provided  f  is  a  solution  with 
i.-O  The  ezpressic:  for  f  is  given  by  Eq.  P.i  as 

U+  1 


r  ,  -  (i  -f)  f 

n+A  'r  dr  r* 


B.3 


i.i  the  Jifferential  equation  becomes  (by  replacing  r*  by  r.+l  in  Eq.  B.2; 


c  ?  (liliii  + 

S  ^  .  2  +  r 

dr 


,  df  , 
n»3  n+i-, 


dr 


)  =  f, 


r»+l  . 


6.4 


labs*.  i*ute  Eq.  B.3  into  Eq-  B  4 

df 


r  .2  ,  df  .  ,  .  df  "I  -.2  ,  df 

2  d  /i _ n  .  2n+3  d  l  n*  _  A /  A_ re 

1-2  r  r  dr  'r  dr  'J  “  -,.2  ^r  dr  ' 

‘•dr  J  dt 
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,  d:'f 

(A - £ 


2  a2fn 


dr" 


_ _2_  2n 

r2  ar2  r’  br)'  r 


a 


dt“~dr 


,2  [l 

'i  U 


a3f 


n  2n+I 


/i£fa  ;La%l  i!\ 

r  dra  r2  dr  -1  r  dt2d: 


,  d2f  of  1 

(£ _ E  _  A  _El 

r  2  2  5iT' 

or  r  J 


Sr*' 


dyf 

_ n 

dt2dr 


B.5 


However 


1 

r 


Sr2 


I_*fn 

2  dr 


s  a  8f 

d_  /I  n-. 
dr  'r  dr' 


Therefore  Eq.  B-5  may  be  written 


r» 


i  d_ 
_r  dr 


a2f 

dr2  * 


.  1  5  ,2n+l  dji.l  _  13  f^E 
r  3r  r  3r  'J  "  r  Sr  ..2 


It  is  apparent  that  this  is  an  identity  since  the  same  expression  may  be 
obtained  by  taking,  tr.e  partial  derivative  of  both  sides  of  Eq.  B.2  with 
respect  to  r  and  dividing  through  by  r.  Therefore  Eq.  B.J  satisfies  the 
differential  equation  and  as  a  consequence  Eq.  B.l  is  a  general  solution  for 
any  value  of  n  provided  f  is  a  solution. 

It  was  shown  in  Chapter  III  by  a  physical  argument  that  f.  has  the 
form  given  by  Eq.  5 .26.  It  is  of  interest  to  show  directly  that  f  in  this 
form  satisfies  the  differential  equation  given  by  Eq  B  2  wi*h  n  -  0.  Tr.;.4- 

is 


+  00 

/  F(t  ‘  et 

o  i 


cosh  u)  du 


is  a  solution  of 


~  d2f  ,  df  d2f 

1  V  *  ^  ’ 


B.6 


B.7 
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provided  that  all  derivatives  which  appear  are  continuous  in  u  and  r.  Take 
the  derivatives  of  Ea.  B.6 


dr  , 

o  -1 

Sr“  ~  C7  . 

oo 

/  F‘  (t  -  g-  cosh  u)  cosh  u  du 

1  ' 

a  i 

a2f 

_ o  _  +1 

dr2  C2 

00 

J'  F"(t  -  g—  cosh  u)  cosh2u  du 
o  1 

32f 

I"”  r 

"T  =  + 

/  F'*(t  -  cosh  u)  du  . 

at2 

o  C1 

Substitute  these  derivatives  into  Eq.  B.7. 

ci[?  J 
C1 

-‘CO 

F"(t  -  7r~  cosh  u)  cosh  u  du 

C1 

0 

«  00 

1 

-  rC, 

J'  F’ (t  -  g—  cosh  u)  cosh  u  duj  =  J'  F"(t  -  g—  cosh  u)  du 

O  1  J  o  1 

Tne  second  tern  on  the  left  may  be  integrated  by  parts  as  follows. 


-  —  cosh  ul  ecsh  u  du  -  +  — ■%—  F' (t  -  cosh  u)  sinh  u 

r°l  vi 


+  -4-  /  F”(t  -  cosh  u)  sinh2u  du 
Jo  l,l 


Tne  first  term,  on  the  right  equals  zero  since  sinh  u  =  0  at  u  =  C  and 

F- (t  -  ecsh  u)  =0  as  u  +  ends  tc  infinity,  for  a  disturbance  which  starts 
C1 


*  In  the  casputat icns  of  Chapters  IV  and  V,  the  function  F"(r|  )  is  not  con¬ 
tinuous  through  7j  -  -a/Cj.  In  this  case  it  is  not  difficult  tc  alter  the 
proof  so  that  the  result  will  hold 
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at  a  finite  time.  Then 


00 

— 7^—  /  F1  (t  -  ~r~  cosh  u)  cosh  u  du 

rCl  Jo  C1 


&  "  “5 


CQ 

^  f *"(*  ~  §- 

C7  J  o  C1 


cosh  u)  sinh  u  du 


O  O 

By  making  this  substitution  into  Eq  B.8  and  observing  that  (cosh  u  -  sinh  u)=  1, 
one  obtains  an  equality. 


00 

'1  [_  "T  S F"^'t  "  C~  cosh  uKcosh2u  ‘  sinh2u)  duj 

03 

=  /  F“(t  -  cosh  u)  du 

Jo  C1 


B.2  General  Form  of  <p  and  V  for  Arbitrary  n. 

The  general  expression  for  the  displacement  potentials  «p  and  as 
given  by  Eq.  5.29)  nay  be  written  in  the  much  more  convenient  form  of  Eq.  5. 30. 
The  derivation  is  shown  for  the  function  cp.  For  V  the  proof  is  the  same. 

From  the  first  of  Eqs.  3  29  >  tne  nth  term  of  the  series  for  <p, 
say  <pn ,  may  be  written 


<Pn  =■  ftn(r.t)  cos  n& 


where 


06 

r  r.,1  A  n  P  r  . 

“n  -  r  (7  J  *'(t-  -  c7  u'  iu’ 


B.  V 


The  term  n  +  1  takes  the  form 

_  n+1-1  d  n+^ 

nJ-1  'r  dr' 


/  F(t  -  5—  cosh  u)  du 
Jo  °1 
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CO 

“n+1  =  r°  I?  (7  IQ  F(t  -  §7  C°Sh  U)  dU‘ 

Tn is  expression  may  be  written  in  the  equivalent  form 


fln+l  -  {  h  Lr”  (7  h]]  -  f  r°  (7  If}  }  7/(t  -  §7  cosh  du- 
By  comparing  this  expression  with  Eq.  B.9  it  may  be  seen  that 


n  .  =  4-  n  -  -  n 

n+1  or  n  r  n 


B.10 


The  derivation  of  Eqs.  3-30  from  Eqs.  3 >29  may  now  proceed  by 

mathematical  induction.  Let  it  be  assumed  that  the  first  of  Eqs.  3. 30  is 

correct  for  0  .  Then  it  will  be  shown  that  Eqs.  3.3O  hold  for  n  , 
n  n+1 

with  n  =  0,  1,  2,  .  .  .  .  Once  this  is  shown  it  only  remains  to  notice 
that  the  expression  is  correct  for  n  =  0.  From  the  first  of  Eqs.  3.30,  q>n 
may  be  written 


<p  -  0  (r  .t)  cos  n&  } 
vn  n'  ’  '  * 


CO 

f l  -  /  Fr\t  -  cosh  u} 

r‘  C?  Jo  C1 


cosh  nu  du  , 


B.  11 


By  performing  the  operation  indicated  by  Eq  B.10,  one  obtains 


jj  -  ~-0 .  /  Fn+-*-(t  -  ~  cosh  u!  cosh  nu  cosh  u  du 

n+!  cn+l  J0  ■  C1 


CO 

—  f  F^t  -  77-  cosh  u)  cosh  nu  du 

r  C?  Jo  C1 


B.12 


The  second  integral  on  the  right  may  be  written  as  fellows  by  integrating  by 


parts 
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cosh  nu  du  = 


f  /  pn(t  ‘  57  cosh  u) 

[r  ‘  5^  COSh  U^n  Sinil  nu^] 


B.13 


OQ 

~  ^AL1-  ^  Fn+1(t  -  cosh  u)(^-  sinh  nu  sinh  u)  du  . 


The  first  term  on  the  right  is  zero,  however,  since  sinh  nu  =  0  at  u  =  0  and 

Fn^,t  •—  cosh  u)  =  0  as  u  tends  to  infinity.  By  using  Eq.  B.13,  one  may 
'“'l 

write  Eq  B.12  in  the  following  form 


(-1) 


n+1 


n+1  n+1 

C1 


f  F“1(t  -  ^ 


cosh  u)  cosh  nu  cosh  u  du 


,  ,\n+l  p  n+1 

«■  -*-= —  /  F  (t  -  ■=—  cosh  u)  sinh  nu  sinh  u  du  , 

C,  Jo  °1 


or  by  combining  terms 


“n.l  -^T -/c  ^  CO*  »)  _ 


ccsh  nu  cosh  u 


t-  sinh  nu  sinh 


uj  du  . 


Tne  q^anti'-y  in  brackets  is,  however,  equal  to  cosh  With  this  sub¬ 

stitution,  the  expression  for  U  .  is  identical  to  Eq.  5-30  with  n+1 

n+± 

replacing  n.  The  conclusion  holds  even  for  n  =  0  since  only  the  first  term  on 
the  right-hand  side  of  Eq.  B.12  appears  in  that  case.  It  has  thus  been  shown 
that  the  first  of  Eqs-  3-30  holds  for  n  +  1  if  it  is  assumed  to  be  correct  for 
n,  n  =■  0,  1,  2,  ...  . 
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To  complete  the  demonstration,  it  is  only  necessary  to  note  that 
the  first  of  E(js.  3-30  is  correct  for  n  =  0,  since  it  is  identical  to  the 
first  of  Eqs.  3'29* 


APPENDIX  C 

DESCRIPTION  OP  THE  COMPUTER  PROGRAM 


A  program  was  written  in  FORTRAN  language  to  solve  the  problem  on 
an  electronic  digital  computer.  The  computer  used  was  a  CDC  1604.  The 
sequence  of  operations  performed  by  the  computer  in  the  solution  of  the 
problem  is  shown  in  the  following  diagrams.  The  computations  are  indicated 
together  with  the  appropriate  equations  to  be  solved. 


Compute  +  radial,  tangential,  uni  hoop  stresses  in 

®the  medium  away  from  the  bcutJLry  usirg  F* "  (rt  ) , 

F"  ' 1  (hi  ) }  .  .  .  .  from  the  above  computations. 
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*  The  superscript  indicates  the  node  (n  -  0,  1,  2).  The  radial  and  shear 
stresses  in  the  diverging  wave  were  also  computed  at  this  time  so  they 
could  later  be  combined  with  the  corresponding  incident  stresses  to 
verify  that  the  sums  were  zero  at  the  boundary. 


Compute  a  new  radius 
r  _  a  +  L(l)^tC1 

a  a 

Compute  the  new  £5,  say  A5‘ 

AS1  =  -  Ab 
r 


Compute  rJ0>  Rjr  4,  .  .  .  . 
\Eqs.  4.2?d,  4.10,  4.12,  4. 26a 
5  =  mA5' ,  m  =  0,  1,  2,  .  .  .  . 


Compute  each  summation  separately  in 
the  expressions  for  stress  given  hy 
Eqs.  4.44  for  n  =  1,  and  similar 
expressions  for  n  =  0  and  2.*  Store 
these  summations. 

P  =  I,  21,  31 . T 


To  obtain  the  stress  in  the  diverging 
wave  combine  the  summations  involving 
F' 1 1  (r^)  with  those  involving  G’  ' 1  (tj^) 
in  each  expression  for  stress  in 
Eqs .  4.44  (for  n  =  l).*  Do  the  same 
for  n  -  0  and  2 . 


Compute  the  incident  stresses. 

Compute  fren  E-i. 

Compute 

(aoV  ^i>p.  '-blV  <ay)p 

....  from 

p  <  2N  +  I (L)  ,  E4S.  4.2-f.f, 

4.23,  4.26f 

p  >  2."  +  I(L) ,  Ecr.  4.25g, 

4.24,  u.2trg 

Combine  outgoing  and  incident  stresses 
for  each  r.cie  ar.i  -r  i*»-  the  total 
stresses  by  cemoining  the  tnree  modes 
for  each  stress  component.  Compute  the 
variation  of  stress  with  angle  around 
the  opening.  Print  these  stresses. 


*  As  discussed  in  the  text,  the  upper  limit  on  the  summations  involving 

G' ' ' (r|2)  actually  used  in  the  machine  was  not  p  -  n/a,  but  p  -  n.  Therefore 
a  shift  in  time  must  be  made  When  these  summations  are  combined  with  those 
involving  F* * ' (up)  in  Eqs.  4.44  and  similar  equations  to  obtain  the  stresses 
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Fig.  j>- j  Vurlm.iun  witn  tliae  of  the  Ft.urit.-r  oorfficicnts  for  stress 
"H  «  oylii.tr -i  1  boundary  during  (■iiV<:jop:r,ent  by  un  incident 
wave  of  '.Ion,  y  3 //,. 
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Fig.  5*3  Comparison  ox’  the  short-time  und  the  computer 
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Fig.  G.l  Hoop  utre;'. a  at  the  boundary  due  to  an  incident  wave  of 
diJutatirn  uoing  n  -  0,  1,  2  (v  «*  0). 
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Fig.  6. j  Comparison  of  the  present,  computations  with  those  of  Re fa .  (3) 
and  (4)  for  hoop  Gtress  at  b  =  0°  due  to  an  incident  wave  of 
dilatation  (v  -  1/5) . 


102 


Time,  tCj/da 

Pig.  6.4  Comparison  of  the  present  computations  with  those  of  Refs.  (3) 
ami  (4)  for  hoop  stress  ut  <>  -*  90°  due  to  an  incident  wave  of 
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